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DEVELOPMENT  OF  CONDUIT  DESIGN 
ANALYTICAL  PROCEDURE 


1 INTRODUCTION 


Background 

Electromagnetic  shielding  is  often  required  in  the  design  of  Army 
facilities  in  order  to  protect  electronic  or  electrical  systems  or 
components  from  the  effects  of  unwanted  electromagnetic  energy.*  These 
effects  range  from  equipment  malfunction  to  actual  component  damage. 
Although  many  sources  of  electromagnetic  interference  (EMI)+  generate 
sufficient  energy  in  stray  or  unwanted  electromagnetic  fields  to  cause 
malfunction  of  electronic  equipment,  damage  from  these  sources  is  un- 
likely. The  possibility  of  damage  is  greatest  from  the  high  energy 
contained  in  an  electromagnetic  pulse  (EMP)  resulting  from  a nuclear 
detonation^  or  the  high  currents  associated  with  a lightning  stroke. 

The  most  common  method  of  providing  protection  against  EMP  is  to 
place  the  equipment  in  electromagnetically  shielded  zones  within  the 
hardened  facility.  In  large,  complex  facilities,  shielding  is  re- 
quired in  several  locations  or  zones,  which  require  electrical  inter- 
connections for  electrical  power  and  data  transmission,  monitor  and 
control  signals,  and  other  general  communication  needs.  The  most 
practical  way  to  make  the  electrical  interconnections  is  by  conven- 
tional wiring  or  cabling.  However,  if  these  wires  and  cables  are 
routed  through  unshielded  zones,  any  EMP  which  occurs  will  induce  po- 
tentially damaging  pulses  on  these  wires. ^ It  is  therefore  essential 


*Two  recent  examples  are  the  SAFEGUARD  BMD  site,  which  required  80-dB 
EMP  shielding  in  addition  to  EMI  shielding  of  certain  portions  of  the 
facility,  and  the  STTF  (Systems  Technology  Test  Facility),  Kwajalein, 
which  required  40-  to  76-dB  EMI  sh’elding. 

f'Electromagnetic  interference"  is  a general  term  cormonly  used  to 
describe  any  type  of  radiated  electromagnetic  energy  which  may  inter- 
fere with  operation  of  electrical  or  electronic  equipment  or  instru- 
mentation. Nuclear  electromagnetic  pulse  is  a specific  type  of  EMI 
having  the  form  of  a single  pulse  of  electromagnetic  energy.  The 
shape  and  magnitude  of  the  pulse  (which  determines  its  spectral  energy 
density  content)  is  predictable  for  various  threats.  Radio  frequency 
interference  is  another  specific  kind  of  EMI  within  the  radio  fre- 
quency spectrum. 

‘B.  D.  Favaudo  and  L.  C.  Martin,  Review  of  Factors  for  Application  in 
Component  Damage  Analysis,  Protection  Engineering  and  Management  (PEM) 
note,  PEM-52  (September  1976). 

^E.  F.  Vance,  Coupling  to  Cables,  DNA  Handbook  Revision,  Chapter  11, 
ADB001204  (Defense  Nuclear  Agency,  December  1974);  and  E.  F.  Vance, 
Design  Guidelines  for  the  Treatment  of  Penetrations  Entering  Communi- 
cations Facilities,  ADB007076  (Defense  Nuclear  Agency,  August  1975). 
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that  all  wiring  and  cabling  penetrating  the  shielded  zones  also  be 
shielded  against  EMP.  Such  shielding  can  be  accomplished  by  using 
individually  shielded  wires  or  cables  or  by  routing  the  wiring  and 
cabling  through  metal  conduits.  Physical  protection  and  shielding  can 
be  provided  by  these  conduits,  which  are  typically  used  for  electrical 
entry^  and  have  been  used  for  signal  cable  routing  in  EMP-protected 
systems . '* 


The  facility  designer  must  therefore  be  aware  of  the  shielding 
properties  of  the  various  components  which  make  up  a conduit  system  so 
that  an  adequate  shield  for  the  wiring  and  cabling  will  be  provided. 

The  components  of  the  conduit  system  may  include  the  conduit,  conduit 
fittings,  and  other  conduit-related  hardware  items  such  as  couplings, 
unions,  junction  boxes,  pull  boxes,  cast  fittings,  and  flexible  joints. 
Most  of  these  related  items  are  standard  electrical  items,  since  few 
items  have  been  developed  specifically  for  EMP-hardened  applications.® 
The  conduit  system  does  not  provide  infinite  shielding,  but  is  com- 
promised by  the  limited  shielding  of  the  conduit  itself,  by  deficien- 
cies inherent  in  some  standard  hardware  items,  and  in  some  cases  by 
improper  installation  practices.® 


To  date,  conduit  design  has  been  accomplished  without  adequate 
guidance  which  would  enable  accurate  determination  of  the  conduit  system 
shielding  and  prediction  of  induced  pulse  energy  on  conduit-housed 
circuits.  The  overall  EMP  protection  acceptability  determination  has 
been  accomplished  after  design  and  construction  by  extensive  analysis 
and  field  and  laboratory  testing.  These  extensive  evaluations  and  other 
studies  in  this  technological  area  have  generated  a wealth  of  technical 
data,  but  using  these  data  requires  the  designer  to  review  an  excessive 
amount  of  literature,  not  all  of  which  is  readily  available.  Although 
TM5-855-5^  contains  a variety  of  EMP  design  information,  it  does  not 
contain  the  detailed  information  required  for  conduit  system  design  and 
evaluation. 


'E.  F.  Vance,  Electromagnetic  - Pulse  Handbook  for  Electric  Fewer 
Systems,  ADA00922R  (Defense  Nuclear  Agency,  February  1975),  p 145. 

‘‘D.  J.  Leverenz,  R.  G.  McCormack,  and  P.  H.  Nielsen,  EMP  Shielding 
Properties  of  Conduit  Systems  and  Related  Harduare,  Technical  Report 
C-1 9/ADAOl 2729  (U.  S.  Army  Construction  Engineering  Research  Labora- 
tory [CERL],  June  1975),  p 11. 

^As-Built  Survey  and  Evaluation  of  EtPP/RFI  Protective  Features,  Vol . 2, 
HND-SP-72-145-ED-R  (Huntsville  Division  [HND]  of  the  U.  S.  Army  Corps 
of  Engineers,  August  1973);  and  SAFEGUARD  EMF’/RFI  Lessons  Learned 
(SAFEGUARD  Ground  Facilities),  HND-SP-75-350-ED-SR  (HND,  31  December 
1975). 

®D.  J.  Leverenz,  R.  G.  McCormack,  and  P.  H.  Nielsen,  Dexielopment  and 
Evaluation  of  Repairs  and  EMP  Leaks  in  Conduit  Systems,  Technical 
Report  C-1 7/ADAOl 1223  (CERL,  April  1975). 

^Nuclear  Electromagnetic  Pulse  Protection,  TM  5-855-5  (Department  of  the 
Army,  February  1974). 
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An  analytical  design  procedure  which  can  be  used  to  theoretically 
predict  the  shielding  of  complex  conduit  systems  is  therefore  needed. 

Such  a procedure  would  allow  the  conduit  system  designer  to  systematical- 
ly assess  the  EMP  hardness  of  conduit  systems  and  assure  a cost-effective 
design  which  provides  adequate  EMP  protection  of  sensitive  electronic/ 
electrical  equipment. 


Objective 


The  objective  of  this  study  is  to  develop  a model  of  conduit  system 
shielding  behavior  for  use  in  an  analytical  procedure  which  is  to  be  in- 
cluded in  TM  5-855-5.  The  procedure,  which  will  be  used  in  the  design 
of  EMP-hardened  conduit  systems,  is  intended  to  provide  District  engineer 
personnel  with  a practical  basis  for  evaluating  the  potential  EMP  pene- 
tration (and  hence  shielding  effectiveness)  of  the  more  common  conduit 
configurations,  and  for  predicting  the  induced  signals  on  conductors 
within  conduit  systems. 


Scope 


This  report  presents  preliminary  models  for  EMP  leakage  into  elec- 
trical conduits.  These  models  are  intended  for  use  in  an  analytical 
procedure  for  conduit  system  design  under  development  at  CERL.  The 
analytical  procedure  will  enable  a designer  to  predict  the  EMP  shield- 
ing performance  on  EMP  hardness  of  proposed  conduit  systems.  The  model 
developed  represents  a significant  portion  of  the  final  procedure; 
however,  for  actual  determination  of  signal  levels  induced  on  conduit 
protected  conductors,  the  following  additional  information  is  necessary: 
(a)  the  current  induced  on  the  conduit,  (b)  cross-coupling  between 
conductors  in  the  same  conduit,  (c)  EMP  leakage  characterization 
(analytical  and/or  empirical)  of  conduit-related  hardware  items.  This 
information  can  either  be  derived  from  the  literature  or  further  ex- 
perimental and  mathematical  studies.  The  final  analytical  procedure 
will  include  enough  of  this  information  to  be  independent  and  not  rely 
on  external  sources  of  information. 

The  model  developed  can  be  applied  to  shielded  cables;  however, 
since  braided  cables  have  additional  coupling  mechanisms  not  present  in 
solid  shields,  the  model  will  not  give  an  accurate  representation  of 
the  shielding  behavior  of  braided  cables. 


Approach 

Development  of  the  leakage  model  for  conduit  systems  involved  three 
major  efforts:  (a)  literature  search,  (b)  mathematical  analysis,  and 
(c)  empirical  studies.  The  detailed  literature  search  was  performed  to 
provide  a summary  of  previous  efforts  directed  toward  development  of 
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mathematical  models  or  procedures  for  predicting  conduit  system  shield- 
ing. Analysis  of  existing  data  then  permitted  determination  of  the 
form  for  general  models  for  predicting  the  signals  on  various  circuits 
within  conduits  for  specific  conduit  defects  and  for  the  inherent  con- 
duit shielding.  Mathematical  analysis  was  performed  to  expand  and  re- 
fine the  models  and  show  expected  conduit  circuit  signals  for  a 
generalized  pulse  threat.  Both  time  domain  and  frequency  domain  ana- 
lyses were  investigated.  Empirical  work  utilizing  pulse  current  in- 
jection testing  was  then  performed  to  verify  the  time  domain  models 
and  determine  values  for  mathematical  model  coefficients  for  the  var- 
ious conduit  hardware  items. 


Organization  of  Report 

Chapter  2 describes  the  model  for  conduit  leakage.  The  appendices 
provide  detailed  information  on  the  development  of  the  model.  Pertinent 
characteristics  of  EMP  are  discussed  in  Appendix  A.  Appendix  B des- 
cribes the  general  experimental  procedures  used  to  obtain  conduit  and 
related  hardware  parameters  necessary  for  the  calculations  used  in  the 
model.  Appendix  C presents  the  theoretical  model  of  EMP  penetration  of 
solid  conduit  and  the  experimental  verifications.  Appendices  D and  E 
describe  typical  values  of  design  dimensions  and  material  properties  of 
conduit,  respectively.  Appendices  F through  H detail  the  theoretical 
analyses  and  experimental  evaluations  used  in  developing  the  models. 

The  annotated  bibliography  presents  references  which,  although  not  re- 
quired for  use  of  the  model,  provide  a source  of  additional  information. 
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DESCRIPTION  OF  THE  CONDUIT  LEAKAGE  MODEL 


Overview 


Efforts  were  directed  toward  developing  and  verifying  mathematical 
techniques  for  determining  the  EMP  field  penetration  into  conduit  sys- 
tems and  the  level  of  induced  signals  on  lines  within  the  conduit  re- 
sulting from  such  penetration.  Although  the  analytical  procedure  occa- 
sionally involves  complicated  mathematical  operations,  computer  calcula- 
tions have  held  to  a minimum.  The  proposed  procedure  uses  simplified 
mathematical  models  which  describe  the  major  features  of  the  induced 
signals,  supplemented  with  empirical  data  from  injected-current  EMP 
tests  of  various  hardware  items.  The  procedure  consists  of  the  follow- 
ing steps: 

a.  Determination  of  conduit  current  and  the  acceptable  signal 
levels  (degree  of  hardness  to  be  maintained).  The  conduit  current, 
including  magnitude  and  time  variation,  must  be  established.  The 
actual  nature  of  the  conduit  current  depends  on  a number  of  factors-- 
often  in  a rather  complicated  manner.  It  is  beyond  the  scope  of  this 
study  to  detail  the  calculation  of  EMP  pickup  by  conduit  systems. 
Appendix  A describes  some  general  characteristics  of  EMP  and  the  bibli- 
ography contains  some  references  on  the  topic  which  can  be  used  by  the 
designer.  Alternatively,  the  threat  current  may  be  supplied  by  the 
basic  construction  specification  or  may  be  assumed  to  have  a peak 
magnitude  of  500  to  10,000  amperes®  with  a waveshape  similar  to  the 
incident  EMP  radiation.  The  degree  of  hardness  required  depends  on  the 
components  and  equipment  connected  to  the  conduit  protected  wires  and 
cables.  This  information  should  be  provided  in  the  installation  speci- 
fications. 

b.  Establishing  the  basic  conduit  design  configuration  and  iden- 
tifying the  type  and  quantity  of  system  components,  including  sections 
of  conduit,  conduit  size,  and  conduit  material;  couplings;  unions; 
unilets,  condulets,  etc;  pull  boxes;  and  junction  boxes.  Any  defects 
such  as  openings  or  electrical  discontinuities  must  also  be  identified. 

c.  Calculation  of  the  EMP  penetration  of  the  conduit  for  the 
specified  current.  Because  the  basic  shielding  is  provided  by  the 
solid  conduit,  the  analyses  start  by  considering  the  EMP  penetration 
through  the  conduit  wall  itself.  The  procedure  for  this  calculation  can 
use  either  time  domain  or  frequency  domain  analysis.  This  calculation 
is  discussed  in  greater  detail  in  the  section  of  this  chapter  entitled 
Calculation  of  EMP  Penetration  of  Solid  Conduit. 


F.  Vance,  Treatment  of  Penetrations  Entering  Communiaations  Facili- 
ties, ADB007076  (Defense  Nuclear  Agency,  August  1975),  p 2. 


d.  Determination  of  the  effects  of  defects  and  components  for 
specified  current.  In  developing  the  mathematical  model  for  the  EMP 
penetration  of  some  defects  and  hardware  items,  emphasis  was  placed 
on  coordinating  simple  mathematical  models  with  results  of  empirical 
evaluations  of  these  items.  This  was  done  because  the  limited  use  of 
some  of  these  items  did  not  warrant  a detailed  mathematical  analysis 
based  on  electromagnetic  field  theory.  Although  testing  every  type  and 
variation  of  every  hardware  item  was  not  possible,  a representative 
selection  has  been  tested  and  their  properties  tabulated.  The  proce- 
dure for  this  calculation  can  use  either  time  domain  or  frequency 
domain  analysis.  This  step  is  discussed  further  in  the  section  of  this 
chapter  entitled  Calculation  of  EMP  Penetration  of  Defects,  Conduit 
Fittings,  and  Related  Hardware. 

e.  Accounting  for  propagation  characteristics  of  the  signal.  For 
rapidly  varying  EMP-induced  signals,  the  propagation  must  be  analyzed 
using  transmission  line  theory.  The  onset  of  the  signal  at  the  end  of 
a conduit  system  will  be  delayed  by  the  time  it  takes  the  signal  to 
travel  the  distance  from  the  location  of  the  defect  to  the  end  of  the 
conduit.  Furthermore,  possible  reflections  due  to  mismatch  in  termina- 
tion of  a transmission  line  may  also  arise. 

For  slowly  varying  signals,  the  transmission  line  properties  be- 
come less  significant,  and  in  the  limit,  such  problems  approach  the 
equivalent  circuit  problem  with  lumped  parameters.  The  time  delay  due 
to  propagation  time  becomes  negligible  for  slowly  varying  signals, 
since  the  velocity  of  propagation  is  nearly  that  of  the  speed  of  light. 

f.  Accounting  for  the  effect  of  circuit  impedances  and  wire  and 
cable  type.  A general  mathematical  model  was  developed  for  expressing 
the  short  circuit  current  in  terms  of  inherent  resistance  (R)  and 
inductance  (L).  This  model  can  be  modified  to  show  the  effects  of  in- 
creasing R and/or  L.  This  step  is  discussed  in  the  section  of  this 
chapter  entitled  Effects  of  Circuit  Configuration  on  EMP-Induced  Signals. 

g.  Accounting  for  possible  signal  summation.  Development  of 
methods  for  performing  this  step  is  planned  for  future  studies. 

h.  Comparison  of  the  composite  EMP  signal  with  the  equipment 
susceptibility. 

i.  Redesign  of  the  conduit  system  as  required,  based  on  problem 
areas  identified  in  step  h. 


Calculation  of  EMP  Penetration  of  Solid  Conduit 

This  step  involves  determining  the  electromagnetic  signal  induced 
in  the  interior  of  a hollow  circular  cylinder  because  of  a current 
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pulse*  on  the  outer  surface.  Since  the  current  pulse  can  take  various 
forms,  all  of  varying  frequency  content,  a general  solution  (assuming 
constant  permeability  and  conductivity)  can  be  determined  by  calculat- 
ing the  time  response  of  the  conduit  due  to  an  impulse  of  electro- 
magnetic energy,  and  then  convoluting  this  response  with  any  given  input. 
Similarly,  the  Fourier  transform  of  the  internal  electric  field  re- 
sulting from  an  impulse  current  (referred  to  in  the  literature  as  the 
surface  transfer  impedance)  can  be  multiplied  by  the  Fourier  transform 
of  the  current  pulse  to  obtain  the  frequency  domain  solution  for  the 
electric  field. 

To  simplify  calculation  of  the  electric  field  inside  the  conduit, 
the  modified  Bessel  functions  appearing  in  the  Laplace  transform  of  t-he 
electric  field  are  replaced  in  this  model  by  their  asymptotic  expan- 
sions. The  results  obtained  using  the  asymtotic  form  of  the  Bessel 
functions  are  valid  for  cylindrical  shields  that  have  inner  radii  which 
are  large  compared  to  the  skin  depth  of  the  lowest  frequency  of  interest 
in  the  surface  current  pulse. + Two  infinite  series  solutions  for  the 
electric  field  due  to  current  impulse  have  been  developed  (see  Appendix 
C)--one  which  converges  for  early  times  (Eq  1)  and  one  which  converges 
for  late  times  (Eq  2): 


E(a,  0)  = Q [• 


ya^(b-a)3  ^ 


«>  exp  [-  ) ] 

z — - 0 ■' 


n=0 


-(2n+l)‘ 


- 1/2] 


e > 0 


[Eq  1] 


*The  term  "pulse"  is  used  to  refer  to  an  arbitrary  current  injection 
of  short  duration;  the  term  "impulse"  is  reserved  for  a signal  of 
negligible  duration  whose  magnitude  is  sufficiently  large  to  induce 
a finite  change  on  the  system. 

■•■It  is  shown  in  Appendix  C that  a (radius)  to  6 (skin  depth)  ratio  is 
10  for  a 1-in.  steel  conduit'at  a frequency  of  200  kHz  assuming 
values  of  150  for  the  relative  permeability  and  20  x 10"°  ohm-m  for 
the  resistivity. 
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E(a,  e)  = Q L — 3/^ 


ya^(b-a)^  ^ 


— — 00  1 2 2 

q [E  (-1)""'  s-i- 

^ n=0 


n2^2 

exp  (-^  0)] 


[Eq  2] 


, 0 > 0 


wnere  E = internal  electric  field 
Q = injected  charge 
a = inner  radius 
b = outer  radius 
a = conductivity 
p = permeability 

0 = relative  time  expressed  in  units  of  a characteristic 
diffusion  time  0 = t/T 

j ^ ya(b-a)2  ' ^ 

' 4 

E(a,  0)  can  be  put  in  the  form 

E(a,  0)  = QFS  (0)  [Eq  3] 

where  Q = charge  passed  along  the  cylinder 

^ ya‘^(  b-a]^ 

0 = t/T 

S (0)  = the  two  series  representations  denoted  by  braces  in 
Eqs  1 and  2. 


Thus,  the  unit  impulse  response  for  a circular  cylinder  of  a 
specified  size  and  material  can  be  defined  by  the  shield  parameters  T 
and  F,  which  depend  on  the  dimensions  and  material  properties  of  the 
cylinder.  Investigation  (Appendix  C)  showed  that  these  two  parameters 
can  be  easily  determined  experimentally,  which  is  important  since  the 
material  properties  of  conduits  are  not  readily  known.  Experimentation 
also  indicated  that  once  the  parameters  for  one  size  cylinder  of  a 
given  material  are  known,  the  parameters  for  another  size  cylinder  of 
the  same  material  can  be  found. 


18 


ssTasaen 


► ' ■ 


i 


I 


i 

i 


The  electric  field  for  a current  pulse  (I  ) other  than  an  impulse 
can  be  determined  by  convolution  with  the  impulse  response  (Eq  1 or 

Eq  2): 

t 

E(a,  t)  = QF  / (t-X)S  (A/T)d  X [Eq  4] 


where  E = internal  electric  field 
Q = injected  charge 

^ M mo^  (b-a? 

= current  pulse 

S = series  representations  of  impulse  response 
t = time 

2 

T = diffusion  time 
X = an  integration  variable 


Appendix  C provides  the  detailed  theoretical  analysis  and  experi- 
mental evaluation  of  the  model.  The  model  assumes  linearity;  however, 
some  deviations  may  occur  at  high  current  levels  (saturation  effects). 
Several  examples  of  the  impulse  response,  as  well  as  results  for  other 
than  impulse  currents,  are  presented  and  compared  with  the  experimental 
results  and  are  in  excellent  agreement. 


Calculation  of  EMP  Penetration  of  Defects, 

Conduit  Fittings,  and  Related  Hardware 

As  previously  indicated,  the  basic  shielding  of  a conduit  system 
is  provided  by  the  solid  conduit  itself.  However,  defects  in  the 
conduit  such  as  cracks  and  breaks  may  seriously  compromise  the  shielding. 
In  addition,  because  the  many  fittings  and  related  hardware  items  (e.g., 
couplings,  unions,  case  access  fittings,  pull  boxes,  junction  boxes, 
flexible  sections)  used  in  conduit  systems  have  not  been  designed  for 
use  in  EMP-hardened  structures,  they  may  also  compromise  the  EMP  shield- 
ing of  the  system.  Thus,  assessing  the  EMP  penetration  of  the  entire 
conduit  system  requires  that  the  contribution  of  various  defects,  con- 
duit fittings,  and  hardware  items  be  evaluated. 

Because  mathematical  analysis  of  even  relatively  simple  defects, 
conduit  fittings,  and  hardware  items  is  a very  difficult  (and  in  some 
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cases  very  nearly  impossible)  task,  some  reliance  must  be  placed  on  em- 
pirical results.  On  the  other  hand,  the  signals  resulting  from  defects 
quite  often  depend  on  the  magnitude  and  waveshape  of  the  applied  current 
as  well  as  the  characteristics  of  the  defect  itself;  thus  complete 
reliance  on  empirical  data  would  necessitate  the  acquisition  of  enormous 
amounts  of  data  to  account  for  the  potential  variations  in  waveshape 
and  characteristics,  of  defects  which  might  occur  in  actual  conduit 
systems.  The  resultant  volume  of  data  would  be  difficult  if  not  im- 
possible to  use,  and  extrapolation  to  situations  for  which  data  were  not 
taken  would  not  be  easy. 

The  approach  used  is  therefore  a compromise  between  the  theoretical 
and  empirical.  It  uses  relatively  simple  mathematical  models  supple- 
mented with  experimental  evaluation  of  certain  parameters.  i 

To  a first  approximation,  the  signals  induced  on  a wire  within  a | 

conduit  due  to  a defect  can  be  considered  to  result  from  the  wire  coup- 
ling with  the  electric  and  magnetic  fields  associated  with  the  current 
in  the  vicinity  of  the  defect.  The  general  model  for  the  direct  field 
penetration  of  a defect  is  proposed  to  be 

= RI  + M ^ [Eq  5] 

where  Vj^  = voltage  on  the  wire 

R and  M = coupling  coefficients  between  the  wire  and  the  " 

electric  and  magnetic  fields,  respectively 

I = current  pulse. 

Appendix  F describes  development  and  testing  of  the  model,  and 
Appendix  G provides  the  coefficients  for  some  selected  defects,  fittings, 
and  related  hardware  determined  through  empirical  testing. 

1 

Determining  the  Effects  of  Circuit 
Configuration  on  EMP-Induced  Signals 

The  signals  induced  on  a wire  inside  a solid  conduit  due  to  an  EMP 
depend  on  the  configuration  of  the  circuit  associated  with  that  wire. 

This  section  describes  two  configurations— open  and  short  circuits. 

Open  Circuit 

The  open  circuit  voltage  induced  on  a wire  within  a conduit  can  be 
expressed  as  the  EMP-induced  electric  field  integrated  along  the  length 
of  the  conduit.  The  relation  for  this  voltage  is 
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After  this  integration,  all  that  was  introduced  to  modify  the  ex- 
pression for  the  internal  electric  field  was  the  length  of  the  conduit. 
Because  this  modification  is  time-independent,  the  Fourier  transform  of 
the  open-circuit  voltage  expression  is  equivalent  to  the  Fourier  trans- 
form of  the  internal  electric  field.  Thus,  the  open  circuit  voltage 
frequency  response  is  equivalent  to  the  induced  electric  field  fre- 
quency response  described  in  the  section  of  this  chapter  entitled 
Calculation  of  EMP  Penetration  of  Solid  Conduit. 


Short  Circuit 

The  short  circuit  current  response  to  an  impulse  current  pulse  can 
be  described  to  a good  approximation  by  considering  a lumped  parameter 
circuit  being  driven  by  the  integral  of  the  electric  field  over  the 
length  of  exposed  circuit.  Eq  7 gives  an  infinite  series  expression 
developed  using  this  model: 


I (t)  = A 


sin  2}/a 


exp(-a0) 


(-1)  n exp( 


-n  TT 


IT 


[Eq  7] 
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where  a = dimensionless  parameter  which  is  the  ratio  of  the 
conduit  diffusion  time  T and  the  L/R  time  constant 

R = inherent  resistance 

L = inductance 

I = length 

e = t/T 

A = H 1 _L 

2Trv^  /f 

Appendix  H describes  the  theoretical  analysis  and  experimental  evalua- 
tion of  these  models,  which  were  found  to  be  in  good  agreement  with 
experimental  results.  Appendix  H also  contains  graphs  which  enable 
rapid  determination  of  peak  current  and  time  to  occurrence  of  the  peak 
current. 


Frequency  Domain  Analysis 

The  preliminary  analytical  procedure  developed  by  CERL  primarily 
uses  time  domain  analysis.  The  analysis  of  EMP  penetration  of  solid 
conduit  and  of  some  conduit  defects  can  also  be  done  in  the  frequency 
domain.  Development  of  transfer  impedances  and  transfer  functions  for 
use  in  frequency  domain  analysis  has  been  done  by  others.  The  results 
of  some  of  this  recent  work*  are  summarized  in  Appendices  C and  F. 


^H.  A.  Roberts,  J.  Capobianco,  and  F.  Aqee,  SAFEGUARD  Buried  Conduit 
Studies  (Harry  Diamond  Laboratories  [HDL],  undated). 


3 CONCLUSIONS 


This  report  has  presented  a preliminary  model  for  the  design  of 
EMP-hardened  conduit  systems  which  Corps  of  Engineers  design  personnel 
can  use  to  evaluate  the  potential  EMP  penetration  of  the  more  common 
conduit  configurations  and  to  predict  the  induced  signals  on  conductors 
within  conduit  systems.  The  model  developed  is  an  essential  portion  of 
an  analytical  procedure  for  conduit  system  design,  but  additional  infor- 
mation is  necessary  for  completion  of  the  procedure. 

This  report  has  also  detailed  the  development  of  the  proposed  model. 
During  the  development,  the  following  conclusions  were  reached: 

a.  The  penetration  of  conduit  by  a current  pulse  can  be  calculated 
directly  from  the  equations  which  govern  the  propagation  of  electro- 
magnetic fields  through  conducting  media.  An  exact  expression  for  the 
Laplace  transform  of  the  electric  field  due  to  an  impulse  current  can  be 
calculated  in  a straightforward  manner  using  standard  Laplace  transform 
theory,  however,  inversion  of  the  exact  transform  is  somewhat  tedious, 
as  it  involves  finding  the  roots  of  an  expression  involving  Bessel  func- 
tions and  their  cross  products.  A good  approximation  which  involves 
asymptotic  expansion  of  the  Bessel  functions  may  be  used  in  many  cases. 

Two  infinite  series  solutions  for  the  electric  field  in  this  case  are 
presented  (one  of  which  converges  for  early  times  whereas  the  other 
converges  for  late  times).  The  electric  field  for  a current  pulse 
other  than  an  impulse  can  be  determined  using  the  convolution  theorem. 

The  convolution  can  be  accomplished  using  numerical  integration.  Sever- 
al examples  of  the  impulse  response,  as  well  as  results  for  other  than 
impulse  currents,  are  presented  and  compared  with  experimental  results 
with  excellent  agreement. 

b.  The  current  arising  on  a shorted  wire  due  to  the  electric  field 
penetration  of  rigid  conduit  was  found  to  be  described  to  a good  approxi- 
mation by  considering  a lumped  parameter  circuit  being  driven  by  the 
integral  of  the  electric  field  over  the  length  of  exposed  conduit.  An 
infinite  series  expression  using  this  model  has  been  developed  and 
found  to  be  in  good  agreement  with  experimental  results. 

c.  The  material  properties  for  construction  grade  conduit  are  not 
readily  available. 

d.  Transfer  impedances  or  transfer  functions  present  an  alternate 
analytic  approach  for  conduit  system  design.  Transfer  functions  for 
solid  conduit  and  a limited  number  of  conduit  or  related  hardware  de- 
fects are  available  from  the  literature  (see  the  annotated  bibliography) 
or  from  work  by  other  agencies.  If  the  frequency  domain  approach  is 
pursued,  measurements  on  additional  items  should  be  made.  Since,  in 
general,  the  transfer  impedance  may  be  a complicated  function,  computerized 
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i 

I numerical  integration  techniques  may  be  necessary  to  obtain  the  system 

j response  in  the  time  domain.  Existing  programs  will  need  to  be  adapted 

[ or  new  programs  developed  for  this  requirement. 

e.  A mathematical  model  can  be  developed  for  conduits  with 
I aperture-type  defects  which  will  define  to  a first  approximation  the 

i induced  current  on  wires  internal  to  the  conduit.  A table  of  coeffi- 

j cients  for  this  model  can  be  prepared  through  empirical  work;  develop- 

i ment  of  such  a table  will  enable  a rapid  general  solution  to  EMP  threat 

1 currents. 
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APPENDIX  A: 

CHARACTERISTICS  OF  EMP 


The  EMP  is  characterized  by  large  magnitude  electromagnetic  fields 
with  short  rise  times  (a  few  nanoseconds)  and  somewhat  longer  fall 
times  (a  few  microseconds).  Most  of  the  energy  associated  with  the  EMP 
lies  within  the  radio- frequency  spectrum  with  a range  from  a few  hertz 
to  the  very  high  frequency  (VHF)  band. 

The  EMP  differs  from  most  other  natural  (e.g.,  lightning)  or  man- 
made (e.g.,  radar)  EMI  in  that  the  EMP's  time  waveform  shows  a faster 
rise  time  and  a higher  amplitude.  Normally,  most  natural  and  man-made 
sources  are  confined  to  a narrow  portion  of  the  frequency  spectrum, 
whereas  the  EMP  occupies  a rather  wide  portion.  In  addition,  high- 
intensity  EMP  fields  can  occur  over  a large  area  almost  simultaneously. 
Intense  natural  and  man-made  electromagnetic  fields  seldom  have  such  a 
wide  simultaneous  distribution. 


EMP  fields  can  have  considerable  variation  of  magnitude,  time 
dependence,  and  direction.  One  can,  however,  define  a typical  field 
environment  which  possesses  the  major  features.  For  the  most  part,  a 
double-exponential  approximation  can  be  used  to  describe  the  most  im- 
portant features 

E(t)  = Eo[e-“t-e-^^]  [-£, 

where  a = the  nulse  decay  time  constant 
6 = the  rise  time  constant 


With  a>>6,  E is  approximately  the  peak  value  of  the  incident  electric 
field  strength.  The  values  of  these  quantities  depend  on  weapon 
characteristics  and  the  location  of  the  observer  relative  to  the  burst 
point;  however,  the  following  values  are  representative  and  can  be 
used  for  purposes  of  illustration; 

Eq  ==  5.2  X 10^  volts/meter 
1 


6.7  X 10"  sec 


3.8  X 10"^  sec 

This  waveform  is  shown  in  Figure  A1 . The  Fourier  transform  of 
this  pulse  is 
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[Eq  A2] 


The  magnitude  of  E(w)  is  plotted  in  Figure  A2.  Thus,  a typical 
EMP  is  characterized  by  a very  fast  rise  time  with  a peak  electric 
field  strength  on  the  order  of  magnitude  of  50  kV/m  and  a duration  of 
approximately  1 psec.  The  pulse  contains  significant  energy  at  fre- 
quencies up  to  100  MHz. 

The  fields  incident  on  the  earth's  surface  from  a high  altitude 
nuclear  explosion  can  be  modeled  as  a plane  wave.  The  magnitude  of  the 
electric  field  strength  E (volts/meter)  in  free  space  is  related  to  the 
magnetic  field  strength  H (amperes/meter)  by  the  impedance  of  free 
space. 


E = OqH  = (377  ohms)  H volts/meter  [Eq  A3] 

Since  the  incident  local  field  is  assumed  to  be  a plane  wave,  the  mag- 
netic field  H has  a time  variation  similar  to  the  electric  field  and 
has  a peak  value 

Hpeak  = “377^  amperes/meter  [Eq  A4] 

The  EMP  will,  of  course,  be  modified  by  propagation  through  a material 
other  than  free  space  (or  air). 

The  exact  electrical  transient  induced  on  a line  or  conduit  system 
by  an  incident  EMP  is  a complex  function  of  location  (buried  or  above- 
ground), length,  terminations,  etc.  This  problem  is  analyzed  in  detail 
in  the  Defense  Nuclear  Agency  EMP  Handbook, and  is  not  covered  in  this 
report.  The  analysis  presented  in  this  report  is  not  dependent  on  a 
particular  waveshape  and  is  applicable  to  any  arbitrary  function.  Thus, 
while  it  is  important  to  know  the  induced  signal  for  any  actual  calcula- 
tion of  conduit  system  response,  it  is  not  necessary  to  know  the  exact 
characteristics  of  this  signal  for  an  explanation  of  the  procedure  for 
the  calculations.  A reasonable  approximation  to  use  for  illustration 
purposes  is  a magnitude  of  0.5  to  10  kA^'  with  a pulse  shape  (rise  time, 
duration)  similar  to  that  of  the  incident  field.  This  general  type  of 
transient  would  be  expected  to  appear  on  a fairly  long  aboveground 
line.  A transient  appearing  on  a buried  line  would  be  expected  to  be 
lower  in  magnitude  with  attenuation  of  the  high  frequency  components, 
resulting  in  a slower  rise-time  pulse. 

’ "^E . F.  Vance,  EMP  (Eleatromagnctio  Pulse)  Handbook,  DNA-2114H-2  (Defense 
Nuclear  Agency,  December  1974),  Chapter  11,  ADB001204.  The  total 
document,  DNA-2114H-2,  is  classified  confidential;  however.  Chapter 
11  is  available  in  an  unclassified  version  from  the  Defense  Documenta- 
tion Center  using  the  referenced  AD  number. 

'*E.  F.  Vance,  Design  Guidelines  for'  the  Treatment  of  Penetrations  Enter- 
ing Cormuniaations  Faeilities , ADB007076  (Defense  Nuclear  Agency, 

August  1975). 
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Figure  AT.  Two-exponential  representation  of  the  high-altitude  EMP 
waveform.  From  E.  F.  Vance,  Eleatromagnetic-Pulse 
Handbook  for  Electric  Power  Systems,  ADA009228  (Defense 
Nuclear  Agency,  February  1975),  p 26. 


O 

UJ 

o 


e> 

< 


E 

a> 

I 

> 

I 

Q 

_J 

U 

Li. 

O 

cr 

i— 

o 

UJ 

_J 

LiJ 


FREQUENCY-MHz 


Figure  A2.  Magnitude  of  the  spectrum  of  the  two-exponential  waveform. 
From  E.  F.  Vance,  Electromagnetic-Pulse  Handbook  for 
Electric  Power  Systems,  ADA009228  (Defense  Nuclear  Agency, 
February  1975),  p 26. 
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APPENDIX  B; 


EXPERIMENTAL  PROCEDURE 


As  part  of  the  development  and  evaluation  of  the  analytical  pro- 
cedure, a program  was  developed  to  experimentally  determine  the  EMP 
penetration  of  conduit  with  and  without  defects  and  of  various  conduit 
hardware  items  using  injected-current  tests.  In  addition,  measuring 
fundamental  material  properties  such  as  permeability  and  conductivity 
was  necessary  because  these  properties  are  not  readily  available  for 
construction-grade  materials. 


Pulse  Injection 


Since  the  principal  concern  of  this  study  is  the  effects  of  elec- 
tromagnetic pulses,  the  experimental  evaluation  was  primarily  concerned 
with  pulsed  current.  Two  pulse-injection  techniques  were  used  to 
study  the  effects  of  EMP: 

a.  Capacitor  discharge,  for  pulses  of  short  duration 

b.  Battery  discharge,  for  pulses  of  long  duration. 

The  measurement  setup  for  each  was  essentially  the  same. 

Measurement  Setup 

Figure  B1  shows  a typical  measurement  configuration  consisting  of 
a pulse  generator  (see  the  following  two  sections),  a test  specimen 
composed  of  a section  of  conduit  and  conduit  hardware  items  to  be  tested, 
a shielded  enclosure  to  protect  the  measurement  equipment  used  to 
measure  the  signals  induced  as  a result  of  the  applied  pulse. 

The  test  specimen  was  installed  as  the  ground  side  of  a parallel 
conduit  transmission  line.  One  end  of  the  test  specimen  was  closed 
with  an  end  cap  to  prevent  EMP  leakage,  and  the  other  was  inserted  into 
a conduit  stub  fastened  to  the  shielded  enclosure.  A sense  wire  of  12- 
gauge  solid  copper  was  fastened  to  the  end  cap,  run  the  entire  length 
of  the  test  specimen,  and  extended  into  the  shielded  enclosure  (see 
Figure  B1 ) . In  some  tests,  each  end  of  the  sense  wire  was  terminated 
in  a resistance  approximately  equal  to  the  characteristic  impedance  of 
the  coaxial  system  consisting  of  the  sense  wire  and  test  conduit. 

The  appropriate  electronic  equipment  was  selected  according  to  the 
signal  to  be  measured.  The  signals  most  often  measured  were  (a)  the 
voltage  across  the  termination  resistor  of  the  sense  wire,  and  (b)  the 
current  flowing  in  the  sense  wire  with  the  sense  wire  shorted;  i.e.. 
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Figure  B1 . Parallel  conduit  transmission  line.  (Test  configuration  for  fast 
risetime  pulse  injection.) 
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zero  resistance  termination.  For  low-level  signals,  an  amplifier  with 
a gain  of  950  was  used.  A Pearson  Model  411  current  probe  was  used  to 
detect  the  short  circuit  current.  The  signals  were  fed  into  an  oscillo- 
scope or  an  ERDAC  III  transient  digitizer  and  recorder. 

The  ERDAC  III  is  a digital  recording  oscilloscope  with  the  capa- 
bility of  capturing  and  storing  one-shot  transients.  It  has  a minimum 
sampling  period  of  1 psec  when  sampling  a single  waveform.  It  also 
has  up  to  four-channel  measurement  capability,  but  the  sampling  period 
per  channel  is  multiplied  by  the  number  of  channels  used.  Thus,  for 
two-channel  usage,  as  was  common  in  deriving  much  of  the  data  of  this 
report,  the  sample  interval  per  channel  is  once  in  1 usee.  Each  sample 
is  stored  with  10-bit  resolution  (1  part  in  1024)  and  the  overall 
amplitude  accuracy  is  0.5  percent  of  full  scale.  The  ERDAC  III  has 
cathode  ray  tube  (CRT)  display  of  the  stored  signal  and  a movable  cursor, 
the  X and  y coordinates  of  which  are  read  out  digitally.  Thus,  any 
point  on  a recorded  curve  may  be  quickly  and  accurately  read  out.  The 
ERDAC  III  also  includes  a tape  recording  feature  which  enables  record- 
ing of  all  stored  transients  on  tape  for  later  analysis  by  computer  or 
other  processor.  Signals  with  rise  times  too  fast  to  be  sampled 
accurately  with  the  ERDAC  III  were  monitored  and  photographed  on  oscillo- 
scopes . 

The  applied  pulse  was  monitored  with  a Pearson  Model  310  current 
probe,  and  the  voltage  across  the  transmission  line  termination  was 
measured  with  a Tektronix  6015  high-voltage  probe. 

Modifications  of  this  basic  measurement  configuration  are  described 
in  the  appropriate  sections  of  this  report. 

Capacitor  Discharge  Current  Injection 

The  source  used  for  fast  rise  time  pulse  injection  was  a high 
voltage  spark  discharge  pulser,  as  shown  in  Figure  B1 . A low  induc- 
tance capacitor  was  charged  up  to  the  spark  gap  discharge  voltage.  When 
the  gap  fired,  the  capacitor  was  discharged  into  the  parallel  conduit 
transmission  line.  The  rise  time  of  the  pulse  produced  was  limited 
mainly  by  the  inductance  of  the  capacitor;  the  pulse  length  was  determined 
by  the  value  of  the  capacitor  and  the  load  resistance  on  the  parallel 
conduit  transmission  line.  The  system  used  at  CERL  was  a repetitive 
free-running  pulser  with  a repetition  rate  determined  by  the  voltage  set- 
ting on  the  high  voltage  power  supply.  The  firing  voltage  was  controlled 
by  adjustment  of  the  spark  gap  and  the  concentration  of  the  sulfur  hexa- 
fluoride (SF,)  in  the  spark  gap.  Typical  pulse  peak  values  were  200  A 
with  5 to  lO^nsec  rise  time  and  approximately  4 psec  e"^  fall  time 
(where  e is  the  base  of  natural  logarithms). 

Battenj  Discharge  Setup 

Slower  rise  time  current  injection  was  accomplished  using  the 
battery  discharge  setup  shown  schematically  in  Figure  B2.  Four 
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automotive-type  lead-acid  storage  batteries  in  series  were  used  as  the 
current  source;  they  provided  a capability  for  up  to  250  A peak  injec- 
tion current.  The  current  was  switched  by  an  automotive-type,  engine- 
starter  solenoid  switch  capable  of  switching  several  hundred  amperes. 

The  solenoid  of  this  switch  was  controlled  by  a separate  12-V  storage 
battery  and  a push-button  switch.  The  push-button  switch  was  electric- 
ally connected  by  means  of  a twisted  pair  of  wires  routed  through  a 
mechanical  feed-through  fitting  to  the  interior  of  the  shielded  en- 
closure, thus  providing  simultaneous  control  of  both  the  current  switch- 
ing and  the  measuring  equipment. 

The  injection  current  was  limited  by  means  of  a series  resistor 
made  from  several  parallel  strands  of  nichrome  wire.  The  thermal  time 
constant  of  the  nichrome  wire  was  measured  and  found  to  be  much  greater 
than  the  time  of  occurrence  of  leakage  and  diffusion  current  phenomena 
within  conduits.  The  room- temperature  resistance  of  the  limiting  resis- 
tor was  0.22  ohms,  thus  establishing  the  peak  current  of  approximately 
225  A.  The  inductance  of  the  loop  consisting  of  the  storage  battery 
bank,  conduit,  limiting  resistor,  and  solenoid  switch  was  approximately 
2.9  pH,  resulting  in  a current  rise  time  of  about  13  psec.  Thus,  the 
applied  current  had  the  form 

-t 

I(t)  = Iq  (1  - e ^o)  [Eq  Bl] 

with 
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Yq  13  Visec 

Experimental  Procedure  - Frequency  Domain 

Harry  Diamond  Laboratories  (HDL)  has  performed  some  frequency 
I domain  measurements.'^  The  laboratory  test  arrangement  used  by  HDL  to 

i investigate  conduit  flaws  (Figure  B3)  was  similar  to  that  used  by  CERL, 

I with  the  following  differences: 

I a.  The  conduit  test  sample  length  was  20  ft  (6.1  m) 

I . b.  The  test  transmission  line  was  a coaxial  line 

f 

I c.  The  sense  wire  was  the  center  conductor  from  an  RG  58 

I coaxial  cable  terminated  in  its  characteristic  impedance 

L '^H.  A.  Roberts,  J.  Capobianco,  and  F.  Agee,  SAFEGUARD  Buried  Conduit 

f Studies  (HDL,  undated). 
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d.  The  shielded  room  to  which  the  conduit  sample  was  attached 
was  a 1-m  cube. 

Flaw  impedances  were  measured  using  a Hewlett  Packard  8000  series  spec- 
trum analyzer.  This  instrument  measured  the  amplitude,  but  not  phase 
of 

15(0))  / = Z^(a))  [Eq  B2] 

where  = the  cable  (or  sense  wire)  current 

I^(a))  = the  conduit  current 

Zp(w)  = the  amplitude  of  the  flaw  impedance. 

The  data  were  plotted  as  a decibel  plot  of  I^(a))  / IpCw)  or  equivalently 
as  20  log  [IZp(aj)/2Z  ],  where  Z is  the  characteristic  impedance  of  the 
sense  wire-conduit  transmission^line.  This  is  because  Zp  is  in  parallel 
with  two  transmission  lines  in  series,  each  with  a characteristic  im- 
pedance of  Z^  (Figure  B4). 

The  measurement  of  the  flaw  impedance  consisted  of  replacing  the 
pulser  in  Figure  B3  with  a signal  source  and  measuring  the  ratio  of 
the  conduit  current  to  the  sense  wire  current  at  a sufficient  number  of 
frequencies  to  determine  the  shape  of  the  flaw  impedance  curve.  HDL 
measurements  with  the  spectrum  analyzer  covered  a frequency  range  from 
near  DC  to  100  MHz.  This  measurement  technique  provided  a dynamic 
range  of  approximately  140  dB,  which  is  sufficient  to  measure  the  flaw 
impedance  or  transfer  function  of  flexible  conduit  (0.015  to  0.035  in. 
[0.381  to  0.899  mm]  wall  thickness),  but  not  enough  to  measure  the 
diffusion  current  transfer  function  for  rigid-walled  conduit. 

Instrumentation  to  measure  amplitude  and  phase  information  (both 
of  which  are  necessary  to  transform  from  the  frequency  domain  to  the 
time  domain)  is  available.  This  information  can  be  continuously  plot- 
ted by  instruments  such  as  the  HP  Model  3575A  Gain/Phase  Meter  (1  Hz  to 
13  MHz)  and  HP  Model  8407A  Analyzer  (100  Hz  to  110  MHz).  A vector 
voltmeter  could  be  used  for  manual  plotting  of  the  amplitude  and  phase 
information. 


Measurement  of  Material  Properties 

In  the  calculation  of  the  penetration  of  conduit  by  EMP  (and  hence 
the  calculation  of  the  basic  shielding  provided  by  the  conduit)  pre- 
sented in  Appendix  C,  it  is  apparent  that  the  specimen's  electrical 
conductivity  a and  magnetic  permeability  y are  two  material  properties 
of  interest.  The  values  of  these  properties  are  tabulated  in  a number 
of  places  for  a variety  of  materials.  The  data  reveal  that  the 


aw  impedance  concept. 


properties  are  strongly  influenced  by  chemical  composition  as  well  as 
heat  treatment.  Unfortunately,  values  of  these  properties  could  not 
be  readily  obtained  from  conduit  manufacturers.  The  importance  of 
these  properties  in  the  calculations  and  the  unavailability  of  the  data 
made  experimental  determination  of  these  properties  for  construction- 
grade  conduit  necessary. 

Since  measuring  a specimen's  resistance  is  often  more  convenient 
than  measuring  its  conductivity,  the  electrical  conductivity  a may  be 
determined  as 


i mho 
p meter 


[Eg  B3] 


where  p = the  resistivity  (ohm-meter) 

The  resistivity  is  determined  by  passing  a direct  current  through 
the  specimen  and  measuring  the  voltage  drop  across  the  specimen.  The 
voltage  V,  due  to  a current  I,  is  given  by 


V = IR 

a 

where  R = the  resistance  of  the  specimen  ~ P 
I = the  length  of  the  specimen 
A = the  cross-sectional  area  of  the  specimen. 

In  the  case  of  conduit,  the  cross  section  is  given  by 


[Eg  B4] 


A = IT  (b^  - a^) 


where  a = 
b = 


the  inside  radius 
the  outer  radius 


Hence, 


V = 


pH 

Tr(b^  - a^ 


[Eg  B5] 


[Eg  B6] 


It  follows  that 


- a^)  } 

I 


[Eg  B7] 


Figure  B5  shows  a schematic  of  the  measurement  technique.  Note 
that  the  points  of  measurement  of  the  voltage  are  between  the  points 
of  current  injection.  This  eliminates  the  measurement  of  voltage  due 
to  contact  resistance. 
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of  the  measurement  setup  for  conductivity. 


APPENDIX  C: 


EMP  PENETRATION  OF  SOLID  CONDUIT 


Approach 

This  appendix  describes  the  development  and  testing  of  a method 
for  determining  the  electromagnetic  signal  induced  in  the  interior  of 
a hollow  circular  cylinder  due  to  a current  pulse*  on  the  outer  sur- 
face. Since  the  current  pulse  can  take  various  forms--an  of  different 
frequency  content--a  general  solution  can  be  determined  by  calculating 
the  time  response  of  the  conduit  due  to  an  impulse  of  electromagnetic 
energy  and  then  convoluting  this  response  with  any  given  input.  Simi- 
larly, the  Fourier  transform  of  the  internal  electric  field  due  to  an 
impulse  current  (referred  to  in  the  literature  as  the  surface  trans- 
fer impedance)  can  be  multiplied  by  the  Fourier  transform  of  the  cur- 
rent pulse  to  obtain  the  frequency  domain  solution  for  the  electric 
field. 

Determination  of  the  electric  field  at  the  inner  surface  of  a 
hollow  circular  cylinder  of  finite  conductivity  subjected  to  an  exter- 
nal surface  current  impulse  involves  the  solution  of  Maxwell's  equa- 
tions for  the  axial  electric  field  inside  the  cylinder  due  to  the 
current  applied  to  the  surface  of  the  conduit.  Straightforward  appli- 
cation of  the  Laplace  transform  technique  leads  to  a transform  of  the 
electric  field  which  involves  modified  Bessel  functions.  Similarly, 
the  Fourier  transform  of  the  electric  field  involves  modified  Bessel 
functions  of  complex  arguments. 

In  the  literature^’  concerning  the  electromagnetic  shielding  ef- 
fectiveness of  cylindrical  structures  under  various  transient  conditions. 


*The  term  "pulse"  is  used  to  refer  to  an  arbitrary  current  injection 
of  short  duration;  the  term  "impulse"  is  reserved  for  a signal  of 
negligible  duration  whose  magnitude  is  sufficiently  large  to  induce 
a finite  change  on  the  system. 

”5.  A.  Schelkunoff,  "The  Electromagnetic  Theory  of  Co-axial  Trans- 
mission Lines  and  Cylindrical  Shields,"  Bell  System  Technical 
Journal,  Vol  XIII  p934),  pp  532-579;  Charles  W.  Harrison,  "Trans- 
ient Electromagnetic  Field  Propagation  Through  Infinite  Sheets," 
IEEE  Transactions  on  Antenna  and  Propagation,  Vol  AP-12,  No.  3 (May 
1964),  pp  319-334;  R.W.P.  King  and  E.  W.  Harrison,  "Cylindrical 
Shields,"  IEEE  Transactions  on  Antenna  and  Propagation,  Vol  AP-9, 

No.  2 (1961),  pp  166-170;  D.  A.  Miller  and  P.  P.  Toulias,  "Penetra- 
tion of  Co-axial  Cables  by  Transient  Fields,"  IEEE  Electromagnetic 
Compatibility  Symposium  Record  (1968),  pp  414-423;  E.  F.  Vance  and 
J.  E.  Nanvicz,  "Internal  Voltages  and  Currents  in  Solid  Shielded 
Cables,"  IEEE  Electromagnetic  Compatibility  Symposium  Record  (1963), 
pp  198-209;  and  J.  L.  Erskine,  "Calculations  of  the  Fields  in  a 
Closed  Cylinder  Resulting  from  an  Electromagnetic  Pulse,"  IEEE 
Electromagnetic  Compatibility  Symposium  Record  (1968),  pp  291-297. 
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the  standard  approaches  lead  to  computational  difficulties,  since  the 
presence  of  the  modified  Bessel  functions  makes  direct  inversion  of  the 
Laplace  transform  to  obtain  a time  domain  solution  difficult.  Evalua- 
tion of  the  Fourier  transform  to  obtain  a solution  in  the  frequency 
domain  is  also  difficult.  Many  of  the  techniques  require  extensive 
numerical  calculations;  however,  a large  class  of  problems  exists  for 
which  accurate  and  rapidly  converging  series  solutions  in  the  time 
domain  (and  a relatively  simple  expression  in  the  frequency  domain)  can 
be  obtained  by  replacing  the  modified  Bessel  functions  in  the  trans- 
form(s)  with  their  asymptotic  expansions  for  large  arguments . ‘ The 
results  obtained  using  the  asymptotic  form  of  the  Bessel  functions  are 
valid  for  cylindrical  shields  that  have  inner  radii  which  are  large 
compared  to  the  skin  depth  of  the  lowest  frequency  of  interest  in  the 
surface  current  pulse.  For  a given  cylinder,  the  impulse  response  so 
calculated  can  then  be  convoluted  with  an  arbitrary  surface  current  to 
obtain  the  resultant  electric  field  on  the  interior  of  the  conduit. 

The  impulse  response  for  a circular  cylinder  of  a specified  size 
and  material  can  be  defined  by  two  shield  parameters  which  depend  on 
the  cylinder's  dimensions  and  material  properties.  These  two  par- 
ameters can  be  easily  determined  experimentally,  which  is  important, 
since  the  properties  of  many  materials  are  not  readily  available.  Once 
the  parameters  for  one  size  cylinder  of  a given  material  are  known,  the 
parameters  for  another  size  cylinder  of  the  same  material  can  be  found. 

The  calculated  impulse  responses  were  checked  experimentally.  The 
application  of  the  theoretical  impulse  response  to  the  determination 
of  the  response  for  a complicated  driving  current  was  demonstrated  by 
predicting  the  response  to  a damped  sinusoid  surface  current  and  com- 
paring it  with  the  measured  signal. 


Theoretical  Analysis 


Time  Domain 

The  calculations  for  the  internal  electric  field  refer  tc  a hol- 
low cylinder  (Figure  Cl)  of  length  inner  radius  a,  outer  radius  b, 
conductivity  o,  permeability  p,  and  permittivity  c.  It  is  assumed  that 
a,  p,  and  c are  isotropic  constants. 

The  current,  Ic(t),  applied  to  the  cylinder  is  assumed  to  be  flow- 
ing in  the  positive  z-direction  and  is  uniform  with  respect  to  z and  ij). 
For  the  purpose  of  the  analysis,  it  is  assumed  that  the  applied  current 
is  an  impulse  of  the  form 

^D.  J.  Leverenz,  W.  Croisant,  and  J.  Verdeyen,  "Electromagnetic  In- 
duced Diffusion  Signals  on  Conduit  Protected  Cables,"  1974  lEEE-EMC 
Symiwsium  Record,  IEEE  CHO  8^'’-7  EMC  (1974). 
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Ic(t)  = Q 6(t)  [Eq  Cl] 


where  Q = the  injected  charge 
6{t)  = a unit  impulse 

Defining  A(r,t)  as  the  magnetic  vector  potential  in  the  z-direction 
and  applying  Maxwell's  equations  for  a conducting  medium  (a  medium  for 
which  cf>>coe)  to  this  model  yields  the  following  differential  equations; 

[Eq  C2] 

E(r.t)  = - [Eq  C3] 

where  the  direction  and  field  comoonents  are  defined  in  Figure  Cl. 
The  initial  condition  A(r,0)  = 0 and  the  boundary  conditions: 


9^A(r,t)  , i 3A(r, 
r 9r 


3A(r,t) 
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[Eq  C4] 


= -B(a.t)  = 0 


3A(b,t)  _ _p/k  _ _ pQ6(t) 

^ 2TTb 


[Eq  C5] 


are  to  be  imposed  upon  Eq  C2. 

Use  of  the  Laplace  transform  for  the  time  variable  along  with  the 
specified  spatial  boundary  conditions  results  in  the  following  expres- 
sion for  the  vector  potential  for  the  region  a < r < b: 

P,  K,(4i^  a)  I (4i^  r)  + a)  K„(/uos'  r) 

A(r,s)  = ^ ! 2 ^ 0 


2iTbv/pas  (/pas  b)  K-j  ( /pas  a)  - K-j  ( /pas  b)  Ii(/p^  a) 


where  I and  K = modified  Bessel  functions  of  the  first  and  second  kind, 
respectively. 

The  transform  of  the  electric  field  is  related  to  the  above  by 


E(r,s)  = - sA(r,s) 


However,  for  r = a,  the  exact  relation  (the  Wronskian)- 


I^(w)K^^l(w)  + I^^^(w)-K^(w)  = - 


[Eq  C7] 


[Eq  C8] 


reduces  the  numerator  of  the  bracketed  expression  in  Eq  C6  to  1/41^  a, 
so  that  at  the  inner  surface  of  the  cylinder  the  transform  of  the 
interior  electric  field  resulting  from  the  current  pulse  on  the  outer 
surface  is 


E(3.s)  - [I^(4i^  b)  K^(/p^  a) 

- ( /p^  b ) I ^ ( /p^  a ) j " ^ 


[Eq  C9] 


Due  to  the  modified  Bessel  functions,  it  is  difficult  to  directly 
invert  the  exact  transform  (Eq  C9);  however,  for  large  arguments,  one 
can  use  the  asymptotic  form  of  the  modified  Bessel  functions 
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+ 


[Eq  CIO] 


I (w)  = exfiM  [1 


4v^-1 

8w 


(4v^-1)(4v^-9) 

2!(8w)^ 


.] 


and 


K (w) 

V 


+ 4v^-1 

8w 


+ (4v^-1)(4/-9l., 
2!(8w)^ 


.]  [Eq  Cll] 


In  view  of  Eqs  CIO  and  Cll,  it  is  apparent  that  for  a sufficiently  large 
argument,  the  series  terms  beyond  the  first  may  be  neglected  in  the 
asymptotic  expansions.  The  next  section  demonstrates  that  this  approxi- 
mation is  valid  provided  the  inner  radius  a is  much  larger  than  the  skin 
depth  for  the  lowest  frequency*  of  interest,  i.e., 

a >>  6 = (2/a)uo)^'^^  [Eq  C12] 


Thus,  the  results  using  this  approximation  apply  to  most  problems-- 
surely  to  the  system  considered  in  the  experimental  section.  (However, 
note  that  this  assumption  does  not  place  a restriction  on  the  ratio 
of  the  wall  thickness  to  the  skin  depth.)  Within  this  framework, 

Eq  C9  simplifies  considerably: 

i(a,s)  = -5-  /H  [Eq  C13] 

TT>/ab  ^ exp[(b-a)/pos]  - exp[-(b-a)/uos] 


or 


E(a,s) 


_9_  /£  

2iT/?b  ° sinh[(b-a)v/uosJ 


[Eq  Cl  4] 


where  the  first  term  of  the  series  Eqs  CIO  and  Cll  have  been  used  in- 
stead of  the  Bessel  functions. 

The  Laplace  transforms  (Eqs  Cl 3 and  C14)  can  be  inverted  by  two 
techniques  to  obtain  two  different  series  representations  of  the  elec- 
tric field: 


*A  more  precise  definition  of  this  statement  is  considered  in  the  next 
section  of  this  appendix. 
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a.  Expand  the  transform  (Eq  C13)  in  a Taylor  series  expansion 
and  invert  term  by  term 

b.  Invert  the  transform  (Eq  C14)  directly  using  the  method  of 
residues  to  obtain  a second  series  solution. 

The  details  of  these  inversion  schemes  are  shown  in  the  annex  along 
with  tabular  data  on  the  time  variation.  For  the  first  scheme,  the 
time  history  of  the  electric  field  is  given  by 
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0 = t/T 

S(9)  = the  two  series  representations  denoted  by  braces  in  Eqs 
C15  and  C16 

A table  of  values  of  S(0)  is  presented  in  the  annex  to  this  appendix, 
and  a plot  is  shown  in  Figure  C2. 

It  is  apparent  that  Eq  Cl 5 converges  rapidly  for  small  values  of 
0,  whereas  Eq  C16  converges  rapidly  for  large  values  of  0,  very  slowly 
for  small  values  of  0,  and  not  at  all  for  0=0.  Eq  C15  is  more  con- 
venient for  generation  of  the  impulse  response,  since  the  first  term 
alone  is  dominant  over  most  of  the  range  of  interest.  For  example, 
at  the  peak  of  S(0),  the  first  term  alone  accounts  for  99.99  percent 
of  the  value  of  the  peak.  The  second  term  in  the  series  contributes 

an  8 percent  correction  to  the  first  at  0 = 1.4,  but,  at  this  rela- 

tive time,  the  field  is  only  10  percent  of  its  peak  value. 

The  peak  of  S(0)  occurs  at  a relative  time  given  by 

9p  = 0.36701  [Eq  C18] 

at  which 

S(6p)  = 0.65604  [Eq  C19] 

Therefore,  the  electric  field  attains  a peak  value 

E(a,0p)  = 0.65604  QF  [Eq  C20] 

at 

t = .36701T  (-21 -| 

These  two  characteristics  of  the  peak  electric  field  are  most  useful 
in  the  experimental  determination  of  the  shield  parameters  T and  F. 

For  an  arbitrary  current  pulse  I(-(t),  the  electric  field  can  be 
determined  by  convolution  with  the  impulse  response.  Hence,  the  elec- 
tric field  for  an  arbitrary  current  pulse  is 


E(a,t)  = QF 


I^.(t-A)  S(A/T)  dX 


[Eq  C22] 
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response  function 


Frequency  Domain 

Although  the  main  objective  of  this  study  is  to  examine  the  time 
variation  of  the  electric  field,  it  is  useful  to  briefly  consider  the 
results  obtained  in  the  frequency  domain  using  the  asymptotic  expan- 
sion of  the  Bessel  functions. 

From  the  definition  of  the  Fourier  transform,  it  is  apparent  that 
the  Fourier  transform  of  the  electric  field  at  the  inner  surface  of 
the  cylinder  can  be  obtained  from  Eq  C9  by  letting  s equal  jw.  Thus, 
at  the  inner  surface  of  a conducting  cylinder  the  Fourier  transform  of 
the  electric  field  due  to  a surface  current  impulse  on  the  outer  sur- 
face is 


" 2^  {Il[(l+j)b/6]K^[(l+j)a/6] 

- Ii[(l+j)a/6]K.,[(l+j)b/6]}''' 


[Eq  C23] 


Evaluating  this  expression  is  difficult  because  of  the  modified 
Bessel  functions  of  complex  arguments;  however,  using  the  first  terms 
of  the  asymptotic  expansions  Eqs  CIO  and  Cll,  the  Fourier  transform  can 
be  expressed  in  terms  of  the  dimensionles  variable,  E: 

C = 8uo  (b-a)^a)/4  = 8(b-a)^/6^  [Eq  C24] 


so  that  the  modulus  is  given  by 

jE(a,u))|  = Q[ ^ ] 

2ttv^  (b-a)a  /^osh/^  - cos/^ 

and  the  phase  by 


[Eq  C25] 


i 

I 


(J)  = tan 


tanh  (A/2)  -*•  tan  (/C/2) 
tanh  {/U2)  - tan  (v^72) 


[Eq  C26]  I 

I 


Figure  C3  shows  plots  of  |E(a,(jj)l  and  ()>  as  functions  of  C-  The 
low-pass  filter  characteristics  of  this  response  are  best  illustrated 
by  examining  the  falloff  at  ^ > 100.  For  instance,  the  relative  re- 
sponse at  ^ = 1000  is  84  dB  below  that  at  ^ = 100. 


To  replace  the  Bessel  functions  I-]  and  K]  by  the  approximations 
I which  retain  only  the  first  term  of  the  asymptotic  expansions,  the 

i 

f 

i 
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RELATIVE  AMPLITUDE  (dB)  PHASE  (degrees) 


RELATIVE  FREQUENCY 


Figure  C3.  Phase  and  relative  amplitude  of  the  Fourier  transform  of 
the  impulse  response  as  functions  of  the  relative  frequen- 
cy (Note  that  the  phase  is  continuously  lagging;  how- 
ever, to  conserve  space  the  curve  is  plotted  between  ± 90° 
since  that  is  the  interval  over  which  tan‘l  is  usually 
defined . ) 
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other  terms  must  be  negligible.  The  other  terms  wi 1 1 be  negligible  if 
/paw  a»l  or,  if  6 is  the  skin  depth  of  the  lowest  freguency  component 
of  interest,  then  a/6  >>  !/»/?.  The  exact  degree  to  which  the  terms 
should  be  negligible  depends  somewhat  on  the  particular  problem  being 
considered.  For  instance,  if  a/6  > 10,  then  the  magnitude  of  the 
second  term  is  less  than  5 percent  of  the  first  term.  Assuming  a rela- 
tive permeability  of  150  and  a resistivity  of  20  x 10"^  ohm-m  for  steel 
conduit,  this  corresponds  to  a frequency  of  200  Hz  for  an  a/6  = 10  for  a 
1-in.  (2.54-cm)  conduit.  However,  note  again  that  the  restriction  is  on 
the  radius  of  the  air-filled  core— not  on  metallic  wall  thickness. 


Experimental  Evaluation 


Time  Domain 


The  experiment  depicted  in  Figure  C4  was  used  to  test  the  utility 
of  the  results  of  the  previous  section.  The  ground  side  of  the  trans- 
mission line  was  the  test  conduit,  which  was  connected  to  a conduit 
stud  welded  into  the  side  of  a shielded  room.  The  pulser  consisted  of 
a free-running  adjustable  spark  gap  and  a low  inductance  capacitor 
mounted  inside  an  aluminum  cylinder  pressurized  with  sulfur  hexafluoride 
(SFg)  gas.  This  pulser  was  capable  of  producing  up  to  35-kV  pulses 
with  less  than  a 10-nsec  rise  time  into  a 100-ohm  load. 


The  test  conduit  was  a standard  10-ft  (3.0-m)  section  and  was  thus 
electrically  short  (nanoseconds)  compared  to  the  decay  time  (micro- 
seconds) of  the  capacitor  being  discharged  into  the  transmission  line. 
Thus,  even  with  a mismatch  between  the  transmission  line  and  the  load 
resistor,  Rq,  the  rise  time  of  If.(t)  was  short  compared  to  its  decay 
time,  which  in  turn  was  short  compared  to  the  time  scales  of  response. 

For  this  reason,  the  conduit  driving  current,  Ic(t),  can  be  expressed 
as 

Ic(t)  = f exp(-t/T)  [Eq  C27] 


where  Q = V^C  = the  charge  on  the  capacitor 

T = RpC  = the  time  constant  of  the  capacitor  discharge 

Vq  = the  breakdown  voltage  of  the  spark  gap 

R„  = the  load  resistor  (conduit  transmission  line  termina- 
° tion) 

In  the  limit,  as  t-^O,  I^(t)  yields 

liS  l^it)  = Q6(t)  C28] 


Figure  C4.  Conduit  test  configuration  with  matched  sense  wire. 


Thus,  with  the  proper  selection  of  t,  the  driving  current  produced  in 
this  experiment  can  be  made  to  approximate  an  impulse  function  and 
provide  a means  of  directly  verifying  the  theoretical  calculations. 

To  measure  the  signals  induced  on  a wire  inside  the  conduit,  a 
sense  wire  was  inserted  into  the  test  conduit  as  shown  in  Figure  C4. 
The  sense  wire  was  connected  to  the  pulser  end  of  the  conduit  through 
a terminating  resistor  (Rp)  with  the  end  of  the  conduit  capped  to  keep 
leakage  signals  from  entering  the  system.  The  sense  wire  ran  the  en- 
tire length  of  the  conduit  and  extended  into  the  shielded  room  where 
it  was  terminated  through  the  load  resistor,  R-] , to  the  grounded  wall 
of  the  shielded  room.  The  voltage  developed  across  R^  was  measured 
using  an  operational  amplifier  (gain,  G)  whose  output  was  displayed  on 
an  oscilloscope.  The  scope  traces  from  which  the  experimental  data 
were  taken  were  photographed. 

To  obtain  the  driving  voltage  for  the  sense  wire,  the  electric 
field  was  integrated  over  the  length  of  the  conduit.  It  was  assumed 
that  the  conduit  was  electrically  short  compared  to  the  dominant  wave- 
lengths of  interest  in  the  internal  electric  field  so  that  the  lumped 
circuit  approximation  in  Figure  C5  could  be  used.  The  experimental 
results  show  that  the  time  scales  of  interest  are  tens  of  microseconds 
for  the  aluminum  conduit,  to  milliseconds  for  the  steel  conduit,  which 
makes  this  approximation  valid  except  for  very  long  cylinders. 


Ez 


R| 

(LOAD  TERMINATION) 


Figure  C5.  Equivalent  circuit  for  calculation  of  conduit  response. 
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When  R2  = 0 or  is  small,  the  voltage  displayed  on  the  oscillo- 
scope will  Be 

Vjt)  = GE{a.t)t  = G QF  5(0)^  [Eq  C29] 

where  G = the  amplifier  gain 
l = the  conduit  length 

Thus  this  voltage  becomes  a direct  measure  of  the  internal  electric 
field. 

Since  S(0)  has  a peak  value  of  0.65604  which  occurs  at  0 = 0.36701, 
a direct  measure  of  the  shield  parameters  F and  T can  be  obtained  from 
the  experiment.  If  the  measured  voltage  has  a peak  value  V which 
occurs  at  a time  tp  from  the  time  of  impulse,  then  ^ 

[Eq  C30] 


[Eq  C31] 


(Since  calculation  of  T requires  only  the  determination  of  the  time  at 
which  the  peak  occurs,  which  can  be  done  quite  accurately  from  scope 
traces,  one  can  easily  expect  a 1 percent  accuracy  for  T. 

Thus,  the  two  shield  parameters  needed  to  use  Eq  Cl 7 can  be  di- 
rectly found  from  experimental  results  without  knowing  the  physical 
or  material  properties  of  the  cylindrical  shield.  As  stated  earlier, 
this  is  important,  since  the  values  of  o and  p are  not  always  known. 
Furthermore,  having  determined  the  shield  parameters  F-]  and  T]  for  one 
size  conduit,  the  parameters  F2  and  T2  for  a second  size  conduit  of 
the  same  material  are: 
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where  a = interior  radius 
b = exterior  radius. 

The  experiments  were  performed  on  standard  aluminum  and  galvanized 
steel  rigid  wall  conduits  with  nominal  diameters  of  1 and  2 in.  (25 
and  51  mm).  Table  Cl  gives  the  values  of  interest  for  the  experimental 
condition  for  each  sample.  Using  the  value  of  tp  from  Table  Cl  and 
Eq  C31 , the  value  of  T can  be  found  for  each  sample  so  that  the  abscissa 
of  Figure  C6  can  be  scaled  to  provide  the  normalized  impulse  response 
function  for  each  sample.  Using  the  values  given  in  Table  Cl,  the 
value  of  the  amplitude  scaling  factor  M can  be  found: 

M = .65604F  [Eq  C34] 


Table  C2  gives  the  experimental  values  of  M and  T obtained  for  the  four 
test  samples.  The  experimental  data  were  normalized  to  a common  peak 
of  1.  Figures  C7  and  C8  show  the  comparisons  for  the  four  test  sam- 
ples and  the  excellent  agreement  which  was  obtained  between  theory  and 
experiment. 

Since  the  values  for  p and  a for  the  aluminum  and  steel  used  in 
the  test  conduits  are  not  known,  directly  comparing  the  experimentally 
determined  values  of  M and  T with  the  theoretical  values  is  not  possi- 
ble. Assuming  that  the  material  properties  are  the  same  for  each  alum- 
inum sample  and  the  same  for  each  steel  sample,  some  comparisons  can 
be  made.  From  the  definition  of  T,  it  is  seen  that: 


T (1-in. 

j25-mmj 

aluminum)  _ T 

T (2-in. 

.51 -mm. 

aluminum)  T 

1-in. 


2-in. 


j25-mm 
51  -mm 


steel 


steel 


[Eq  C35] 


which  is  the  case  for  the  values  presented  in  Table  C2.  Based  on  the 
conduits'  dimensions,  from  Table  C3  this  ratio  should  be  1.32,  which 
compares  closely  with  the  experimentally  determined  value  of  1.43. 
This  is  well  within  the  accuracy  which  can  be  expected  due  to  varia- 
tions in  thickness  and  material  properties. 


In  order  to  test  Eq  C22  for  accuracy  in  predicting  conduit  response 
for  a non-impulse  driving  function,  a 1-in.  (25-mm)  aluminum  test  conduit 
was  driven  with  an  exponentially  damped  sinusoid  current.  The  measured 
current  on  the  conduit  was  of  the  form 


I(t)  = exp(-  t/r)  sin(27rft) 

_5 

where  t = 4.62  x 10  sec 

f = 31.5  kHz 


[Eq  C36] 
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Experimental  Constants  for  the  Calculation  of  T and  M for  the  Four  Test  Conduits 


Table  C2 


Experimentally  Determined  Values  of  T and  M 
for  the  Four  Test  Samples 


Test  Sample 

T 

M 

1-in. (25-mm) 
aluminum  conduit 

0.084  msec 

2.546 

2-in. (51-mm) 
aluminum  conduit 

0. 1 20  msec 

0.824 

1-in. (25-mm) 

galvanized  steel  conduit 

4.200  msec 

0.175 

2-in. (51-mm) 

galvanized  steel  conduit 

6.00  msec 

0.0472 
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TIME  IN  MILLISECONDS 


MILLISECONDS 


Table  C3 


I 


I 

i 


t 


Using  Eq  E22  and  the  experimental  value  of  T from  Table  C2,  the  nor- 
malized voltage  across  R]  can  be  predicted  and  compared  with  the  mea- 
sured signal.  The  comparison,  shown  in  Figure  C9,  demonstrates  the  close 
agreement  between  the  predicted  and  measured  responses. 

Ft^equenay  Domain 

Flaw  impedance  or  the  transfer  impedance  of  conduit  is  near  the 
maximum  dynamic  range  of  measurement  instrumentation.  HDL  did,  how- 
ever, make  such  a measurement  on  two  sections  of  flexible  conduit  with 
wall  thicknesses  of  0.015  and  0.035  in.  (0.381  and  0.889  mm).  The 
length  of  the  flexible  section  was  approximately  20  in.  (0.5  m).  Figure 
CIO  shows  the  results  of  this  measurement.  The  theoretical  limit  for 
the  maximum  of  the  transfer  function  is  the  ratio  of  the  DC  resistances 
of  the  conduit  and  the  conductor.  The  flaw  impedance  shows  that  only 
the  low  frequency  components  of  a current  exterior  to  the  conduit  will 
appear  on  a conductor  inside  the  conduit. 


Summary 

It  has  been  shown  that  calculation  of  the  electric  field  on  the 
inner  surface  of  a hollow  circular  cylinder  due  to  a current  pulse  on 
the  outer  surface  can  be  greatly  simplified  with  little  loss  in  accur- 
acy by  replacing  the  modified  Bessel  functions  appearing  in  the  Laplace 
transform  of  the  electric  field  by  their  asymptotic  expansions  for 
large  arguments.  The  results  are  valid  for  those  cylinders  whose  inner 
radius  is  large  compared  to  the  skin  depth  of  the  lowest  frequency 
component  in  the  surface  current  pulse. 

The  calculated  electric  field  due  to  an  impulse  current  has  rela- 
tively simple  series  solutions  which  can  be  used  in  a convolution  pro- 
cess for  an  arbitrary  forcing  function.  It  was  shown  that  these 
electric  field  calculations  can  be  uniquely  specified  by  two  parameters 
which  depend  only  on  the  physical  dimensions  and  material  properties 
of  the  cylinder.  It  was  demonstrated  experimentally  that  these  two 
parameters  are  easily  measured  and  that  experimental  data  can  be  put  in 
the  form  of  the  theoretical  calculations  with  excellent  agreement  for 
1-  and  2-in.  (25-  and  51-mm)  diameter  aluminum  and  steel  cylinders. 
Finally,  it  was  verified  experimentally  that  the  electric  field  could 
be  predicted  for  a complicated  current  pulse  through  convolution  with 
the  impulse  response. 
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1-in.  (25-tiifii)  aluminum  conduit 
on  current. 


FLAW  IMPEDANCE  THIN  WALL  FLEX- JOINT 


Flaw  impedance  of  two  flex-joints  of  different  wall  thickness  without  copper  straps 
(From  H.  A.  Roberts,  J.  Capobianco,  and  F.  Agee,  SAFEGl  Buried  Conduit  Studies 
[HDL,  undated]). 


ANNEX  TO  APPENDIX  C: 

INVERSION  OF  THE  LAPLACE  TRANSFORMS 


The  Laplace  transforms  (Eqs  Cl 3 and  C14)  can  be  inverted  by  either 
of  two  techniques  to  obtain  series  representations  of  the  electric  field: 

a.  Expand  the  transform  (Eq  C13)  in  a Taylor  series  expansion  and 
invert  term  by  term  to  obtain  the  series  expression  of  Eq  CIS. 

b.  Invert  the  transform  (Eq  C14)  directly  using  the  method  of 
residues  to  obtain  a second  series  solution  given  by  Eq  C16. 

To  apply  the  first  method,  Eq  Cl 3 is  rewritten 


p/,  s)  = r Q /H-i  exp[-(b-a)  /yos] 
fi/Jb  1 - exp[-2(b-a)/uosJ 


[Eq  C37] 


using  the  Taylor  series  expansion, 


= 1 + X + X + 


= I x" 

n=0 


[Eq  C38] 


Then; 


E{a,s)  = [-^ — /£]  /s  exp[-(2n+l ) (b-a)/pasj 
TTv/aF  ^ n=0 


[Eq  C39] 


The  inverse  for  each  term  is 


, 1 exp[-kj^/4t]^k^  1^ 


[Eq  C40] 


If  a parameter,  T,  is  defined  as  depending  only  on 

material  properties  p and  a and  the  wall  thickness  (b-a),  then  a 
dimensionless  variable,  0,  may  be  defined  as  t/T  so  that  the  electric 
field  at  the  inner  surface  of  a hollow  cylinder  caused  by  an  impulse 
current  on  its  outer  surface  may  be  expressed  as 


E{a,e)  =Q[- 


3^- 

■/ab  pa  (b-a) n=0 


exp[-! 


‘[12^.^])  PEq  C41] 


63 


i 


y j 

t Ur*’ 


Table  C4  evaluates  the  term  in  braces  in  Eq  C41 . 

To  obtain  the  second  series  representation  consider  Eq  C37  in  the 


form 


E(a,s)  • [-L  /f]  - 

n _ /-» u y o 


2TT/ab  ^ sinh(2/r  /s) 


[Eq  C42] 


then  the  inverse  is 

E(a,s)  = [ 


_JL_  /£]  _J_ 
/H  “ 


exp(st) 


Y-J 


. sinh{2/r  /s) 


ds 


[Eq  C43] 


= [^^  y^]  I R 

27rv^  ° n=0  " 


where  = the  residues  at  the  poles  (sn)  of  the  integrand 
The  residues  are  given  by 


n ^ 

ds 


7s  exp ( St) 


[sinh(2/f  /s] 


s=s. 


_L  A exp(st) 

/f  cosh(2/r  i^) 


[Eq  C44] 


s-s. 


where  s„  is  the  nth  zero  of 
n 


sinh  (2/f~ 7 s ) 


[Eq  C45] 


Since  the  roots  of  Eq  C45  are 


= 


2 2 
n IT 

"4T~ 


[Eq  C46] 


and 


cosh(jmT)  = cos(mT)  = (-1) 


[Eq  C47] 
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Evaluation  of  S(0) 


Table  C4 


exp[- 


.3/2 


,-120+1)5 

L Q 


0/n 

0 

1 

2 

s(e) 

.0 

.00000 

.00000 

,1 

.01364 

.01364 

.2 

.33900 

.33900 

.3 

.61513 

.61513 

• .36701 

« • 

.65604 

.65604 

.4 

.64894 

.64894 

.5 

.57418 

.00000 

.57418 

.6 

.47413 

.00001 

.47414 

.7 

.37997 

.00006 

. 38002 

.8 

.30030 

.00020 

.30050 

.9 

.23562 

.00051 

.23612 

1.0 

.18394 

.00105 

.18499 

1.1 

.14286 

.00186 

.14472 

1.2 

.11020 

.00295 

.11315 

1.3 

.08417 

.00427 

.08843 

1.4 

.06333 

.00578 

.06911 

1.5 

.04658 

.00742 

.05400 

1.6 

.03306 

.00913 

.04219 

1.7 

.02211 

.01086 

. 00000 

.03297 

1.8 

.01320 

.01256 

.00001 

.02576 

1.9 

.00594 

.01418 

.00001 

.02013 

2.0 

.00000 

.01571 

.00002 

.01573 

t 
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then 


, 22  22 
R - H)  n TT  r n TT  I 
\ 4j  expL  tj 


[Eq  C48] 


Thus,  the  second  series  representation  of  the  electric  field  is 


E(a,e)  = Q[- 


8 


ya^(b-a)^ 


] 


7“  “ 1 2 2 

■i'TTr  V f ^ Nfl+I  on  r on 


ij-iy 


2 2 


[Eq  C49] 


n=l 


exp[- 


0]} 


Table  C5  evaluates  the  quantity  in  braces  in  Eq  C49. 


Erskine'®  applied  this  technique  directly  to  the  expression  in 
brackets  in  Eq  C9.  This  latter  expression  is  similar  to  that  obtained 
by  Erskine,  but  without  the  computational  difficulties  of  evaluatinq 
the  residues  the  inversion  integral  of  Eq  C9.  For  thin  wall  cylin- 
ders (b/a=l)  Erskine's  results  reduce  to  an  expression  similar  to  Eq  C16. 


'=’J.  L.  Erskine,  ''Calculation  of  the  Fields  in  a Closed  Cylinder  Result- 
ing from  an  Electromagnetic  Pulse,"  IEEE  Eleatvomxgnetia  Compatibility 
Symposium  Record  (1968),  pp  291-297. 
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Table  C5 


* 


1 


APPENDIX  D: 

DESIGN  DIMENSIONS  AND  PARAMETERS  OF  ELECTRICAL  CONDUIT 


The  EMP  penetration  of  solid  conduit  has  been  shown  to  depend  upon 
the  physical  dimensions  of  conduit  (Appendix  C).  This  appendix  sum- 
marizes the  design  dimensions  for  electrical  conduit.  Table  D1  gives 
the  design  dimensions  for  rigid  metal  conduit,  Table  D2  covers  inter- 
mediate metal  conduit,  and  Table  D3  covers  electrical  metallic  tubing 
(EMT).  Table  D4  summarizes  the  nominal  dimensions  and  several  derived 
parameters  of  ferrous  metal  and  aluminum  that  enter  into  EMP  calcula- 
tions.* 


*Note  that  since  electric  field  is  usually  expressed  in  volts/meter, 
calculations  use  SI  units. 


ii. 
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Table  D1 


Table  02 


Nominal  Dimensions  and  Parameters  of  Ferrous  Metal  and  Aluminum  Rigid  Conduit 


APPENDIX  E: 


MATERIAL  PROPERTIES  OF  ELECTRICAL  CONDUIT 


The  material  properties  of  electrical  conduit  significantly  affect 
the  EMP  penetration  of  the  conduit.  To  make  accurate  EMP  penetration 
predictions,  it  is  necessary  to  know  the  electrical  conductivity,  o, 
and  the  magnetic  permeability,  p. 

Most  electrical  conduit  standards  (and  thus  most  conduit  manufac- 
turers) are  concerned  primarily  with  mechanical  properties  (e.g., 
strength,  cutting  properties,  and  bending  properties)  and  corrosion 
resistance,  so  that  in  many  cases  it  is  difficult  to  find  values  or  to 
deduce  a range  of  values  of  the  electrical  properties  for  electrical 
conduit. 

Since  electrical  conduit  is  not  usually  manufacturered  from  con- 
ductor grade  materials,  the  electrical  conductivity  is  rarely  measured 
or  specified.  Similarly,  the  permeability  is  not  available  for  elec- 
trical conduit  because  the  materials  used  in  conduit  are  not  generally 
regarded  as  electromagnetic  shielding  materials.  However,  although  the 
physical  properties  of  electrical  conduit  are  not  tabulated  in  the 
general  literature,  the  values  of  the  physical  properties  for  closely 
related  materials  may  in  some  cases  permit  an  estimate  of  the  conduit 
values.  This  approach  is  hampered  somewhat  by  variations  in  chemical 
composition,  heat  treatment,  and  mechanical  working--all  of  which  af- 
fect the  electrical  properties. 

Guidance  for  installation  of  electrical  metallic  tubing,  rigid 
metal  conduit,  and  intermediate  metal  conduit  is  provided  by  the 
National  Electrical  Code  (NEC);^®  however,  no  soeci f ication  is  given 
for  conduit  material.  Specifications  for  conduit  manufacture  have  been 
set  by  other  agencies  such  as  American  National  Standards  Institute 
(ANSI)  and  Underwriters'  Laboratories  (UL). 

To  be  approved  by  UL,  individual  conduit  product  lines  (i.e., 
conduit  fittings,  and  junction  boxes)  must  be  manufactured  in  accor- 
dance with  UL  standards.  The  requirements  for  electrical  metallic  tub- 
ing stated  in  UL  797'^  are: 

2.1  Each  tube  used  in  the  manufacture  of  electrical  metallic 
tubing  shall  be  of  mild  steel,  a silicon-bronze  alloy  containing 
at  least  1.25  percent  silicon,  an  aluminum-base  alloy  containing 
no  more  than  0.40  percent  copper,  or  other  metal  which,  upon 

^"^citTohaV  EjectrfcaTl  Code,  1975  Edition,  NFPA  No.  70-1975,  National 
Fire  Protection  Association  (1974).  This  has  been  approved  by  the 
American  National  Standards  Institute  and  is  also  know  as  ANSI 
Cl-1975. 

^"^Standard  for  Eleotriaal  Metallic  Tubing,  UL797  (UL,  Inc.,  June  20, 
1973).  Approved  as  ANSI  33.98-1973,  June  20,  1973. 
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investigation,  is  found  to  be  suitable  for  the  purpose.  The 
tube  shall  be  uniformly  thick  throughout  the  length  of  the  tube. 
The  welding  of  all  seams  shall  be  thoroughly  well  done. 


r 


i 


I 


The  requirements  for  rigid  metal  conduit  stated  in  UL  6^*  are: 

2.1  Each  tube  used  in  the  manufacture  of  rigid  metal  conduit 

shall  be  of  steel,  wrought  iron,  silicon-bronze  alloy  C65100 
(low-silicon  bronze  B)  or  C65500  (high-silicon  bronze  A)  as 
detailed  in  the  American  Society  for  Testing  and  Materials 
"Standard  Specification  for  Copper-Silicon  Alloy  Seamless  Pipe 
and  Tube"  (ASTM  B315-1976),  and  aluminum  alloy  containing  no 
more  than  0.40  percent  copper,  or  another  metal  that  upon  in- 
vestigation, is  found  suitable  for  the  purpose.  ‘ 

The  requirements  for  rigid  aluminum  conduit  specified  in  WW-C-00540c" 
are : 

3.2  Detail.  Conduit,  couplings,  elbows,  and  nipples  furnished 
under  this  specification  shall  conform  to  the  United  States  of 
America  Standards  Institute  (USASI)  Standard  C80.5. 

However,  concerning  the  composition  of  rigid  aluminum  conduit,  USAS 
C80. 5-1966  states:^” 

3.4  Alloy.  The  conduit  shall  be  made  of  an  aluminum  alloy  con-  ' h 

taining  not  more  than  0.40  percent  copper. 

Except  for  the  0.40  percent  copoer,  the  alloy  is  not  specified;  however, 

the  Aluminum  Association  states:^*  I 


Composition  and  Manufacture.  Rigid  aluminum  conduit  is  extruded 
from  the  magnesium-silicide,  6063-Tl  alloy.  Special  precautions 
are  taken  to  obtain  a lower  copper  content  than  industry  stan- 
dards permit.  In  this  way  the  corrosion  resistance  of  the  alloy 
is  greatly  increased. 

Typical  composition  limits  for  aluminum  conduit  metal  are  given  in 
Table  El. 

Nominal  values  for  the  electrical  resistivity  are  shown  in  Table  E2 
for  various  tempers.  This  suggests  that  even  with  the  same  nominal 
chemical  composition,  the  electrical  resistivity  may  vary  up  to  13  per- 
cent depending  upon  the  temper. 


TB 

1 f 


2 0 

2 1 


Standard  for  Rigid  Metal  Conduit,  UL6,  8th  Ed.  (UL,  1°76),  p 5. 

Interim  Federal  Specification  Conduit,  Metal,  Rigid',  and  Coupling , 

Elbow,  and  Ripples,  Electrical  Conduit:  Aluminum,  WW-C-00540c(GSA-FSS) 
(Federal  Supply  Service,  General  Services  Administration,  1967),  p 2. 

USA  Standard  Specification  for  Rigid  Aluminum  Conduit,  USAS  C80.5- 
1966  (ANSIO). 

Aluminum  Electrical  Conductor  Handbook  (The  Aluminum  Association),  p 17-2. 
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Table  El 


Composition  Limits  of  Aluminum  Conduit  Alloy 


Percentage  Limits 
per  Industry 
Standards 

Aluminum 
Conduit,  Typical 
Percentages 

Copper 

0.10  max.* 

0.02 

Silicon 

0.20  to  0.6 

0.40 

Iron 

0.35  max. 

0.20 

Magnesium 

0.45  to  0.9 

0.7 

Manganese 

0.10  max. 

0.01 

Chromium 

0.10  max. 

0.01 

Titanium 

0.10  max. 

0.01 

Zinc 

0.10  max. 

0.02 

Others 

0.15  max. 

Trace 

Aluminum 

Remainder 

Remainder 

♦Maximum  limit  set  by  the  Aluminum  Association.  Alloys  with  up  to  0.40 
percent  copper  are  acceptable  to  UL  for  use  in  rigid  aluminum  conduit; 
typical  conduit  uses  only  0.02,  or  95  percent  less.  (From  Aluminum 
Eleatriaal  Conductors  Handbook,  p 17-6.) 
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Table  E2 


Variation  of  Typical  Values  of  Electrical  Resistivity 
of  Alloy  6063  with  Temper* 


Temper 

Electrical  Resistivity($>m) 

Electrical 

Conductivity  (mho/m) 

0 

3.0X10"^ 

3.4X10’“ 

3.1X10’“ 

3.3X10’“ 

3.3Xlo! 

T1 

2.9X1o4 

T5 

3.2XIO7 

T6,  T83 

3.0x10' 

*From  Aluminum  Standards  and  Data,  5th  Ed.  (Aluminum  Association,  1976), 
p 39. 
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Table  E3  shows  nominal  values  of  resistivity  for  pure  aluminum  cited 
in  the  literature.  Values  of  resistivity  for  aluminum  are  occasionally 
suggested  as  suitable  for  use  in  shielding  studies;  however,  it  is  appar- 
ent that  the  typical  values  cited  for  pure  aluminum  differ  from  those 
given  for  aluminum  conduit  alloy. 

As  part  of  this  study,  the  resistivity  of  one  specimen  of  nominal 
1-in.  (25-mm)  rigid  aluminum  conduit  was  measured.  It  was  found  that 


P 


= 3.8  X 10"® 


Jl-m 


[Eg  El] 


or 


a = 2.6  X 10^ 


mho/m 


[Eg  E2] 


Since  the  alloy  in  aluminum  conduit  is  primarily  aluminum,  one 
would  expect  that  the  permeability  should  be  near  that  of  pure  aluminum, 
i.e.,  a relative  permeability  which  is  essentially  constant  and  nearly 
egual  to  1.  Aluminum  is  actually  paramagnetic  with  a maximum  relative 
permeability^^  of  1.00000065.  The  actual  permeability  of  aluminum 
conduit  alloy  does  not  appear  to  be  available  in  the  literature  or  in 
any  of  the  standards. 

However,  it  is  easier  to  establish  parameters  for  aluminum  conduit 
than  for  ferrous  conduit  because  the  material  composition  is  primarily 
aluminum  with  only  small  guantities  of  alloying  elements  and  impurities. 
The  specifications  for  rigid  steel  conduit,  zinc  coated,  are  given  in 
UL6,  WW-C-581d,^^  and  ANSI  C80. 1-1971  however,  the  composition  of 
the  steel  is  not  specified.  The  reguirements  for  electrical  metallic 
tubing  (EMT)  are  given  in  UL  797,  WW-C-563A,^^  and  USAS  C80. 3-1 966;^® 

Reference  Data  for  Radio  Engineers,  6th  Ed.  (H.  W.  Sams  & Co.,  1975), 
p 4-32. 

^^Federal  Specification  Conduit,  Metal,  Rigid;  and  Cohopling,  Elbow,  and 
Nipple,  Electrical  Conduit:  Zinc  Coated,  WW-C-Bld  (Federal  Supply 
Service,  General  Services  Administration,  5 June  1962). 

^''American  National  Standard  Specification  for  Rigid  Steel  Conduit,  Zinc 
Coated,  ANSI  C30. 1-1971  (R-1966)  (ANSI). 

^^Federal  Specification  Conduit,  Metal,  Rigid:  Electrical,  Thin-Wall 
Steel  Type  (Electrical  Metallic  Tubing);  Straight  Lengths,  Elbows,  and 
Bends,  WW-C-563A  (General  Services  Administration,  December  1973). 
^^USA  Standard  Specification  for  Electrical  Metallic  Tubing,  Zinc 
Coated,  USAS  C80.3-1966  (ANSI). 
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Typical  Values  of  Electrical  Resistivity 
and  Conductivity  of  Aluminum 


p(ohm-meter) 

a(mho/meter) 

Reference 

2.50  X 

10‘®  (0°C) 

4.00  X 

10^ 

★ 

2.74  X 

10"®  (22°C) 

3.65  X 

10^ 

★ 

2.7809 

X 10'®  (20°C) 

3.5960 

X 10^ 

★ * ** 

2.62  X 

10"®  (20°C) 

3.82  X 

10^ 

t 

2.65  X 

10-8 

3.77  X 

10^ 

tt 

*D.  E.  Gray,  Ameriaccn  Institute  of  Physios  Handbook,  3rd  Ed.  (McGraw 
Hill , 1972),  p 9-39. 

**E.  G.  Fink,  Editor-in-Chief , and  J.  M.  Carroll,  Associate  Editor, 
Standard  Handbook  for  Eleotriaal  Engineers,  10th  Ed.  (McGraw  Hill, 
1968),  p 4-8. 

^Standard  Handbook  for  Eleatrioal  Engineers , Reference  Data  for  Radio 
Engineers,  6th  Ed.  (H.  W.  Sams  & Co.,  1975),  p 4-21. 
ttT.  Lyman,  Ed.,  Metals  Handbook,  Vol  1,  8th  Ed.  (American  Society  for 
Metals,  1961 ). 


however,  again  the  composition  of  the  steel  is  not  specified.  In  both 
cases,  requirements  are  made  for  mechanical  properties,  but  the  heat 
treatment  or  mechanical  working  of  the  materials  is  not  specified.  It 
would  be  expected  that  there  will  be  some  variation  in  the  chemical 
composition,  mechanical  working,  and  heat  treatment  of  different  con- 
duit. 

Conduit  conductivity  and  permeability  are  rarely  measured.  How- 
ever, conduit  and  line  pipe  steels  are  generally  similar,  and  the  con- 
ductivity of  line  pipe  steel  has  been  determined.^’  The  resistivities 
(in  micro-ohm-centimeters)  of  24  line  pipe  steels  reoresentati ve  nf  the 
range  of  composition  used  in  manufacturing  API  5L,  API  5LX  (through 
X60),  ASTM  A 53,  and  ASTM  A 106  grades  were  determined.  The  resistivi- 
ties range  from  14.97  to  22.83  uohm-cm. 

Table  E4  shows  the  chemical  compositions  of  the  line  pipe  steel 
specimens.  The  range  (in  percent  by  weight)  of  each  constituent  was: 


Carbon 

Manganese 

Phosphorus 

Sulfur 

Silicon 

Vanadium 


0.13-0.29 

0.39-1 .28 

0.007-0.073 

0.018-0.034 

0.026-0.250 

0.0-0.08 


A multiple  regression  analysis  indicated  that  97.2  percent  of  the  varia- 
tion in  resistivity  could  be  accounted  for  by  the  elements  manganese, 
silicon,  and  phosphorus,  whereas  the  effect  of  carbon,  vanadium,  and 
sulfur  (within  the  range  of  the  steels  considered)  was  not  significant. 
The  resistivity  (in  micro-ohm  centimeters)  can  be  calculated  from  the 
chemical  composition  by  Eq  E3  using  the  percents  of  manganese,  silicon, 
and  phosphorus: 

p = 11.59  + 5.43  (%Mn)  + 16.1  (%Si)  + 15.4  (%P) 


Calculated  and  actual  resistivities  are  compared  in  Table  E5. 


The  permeability  of  steel  is  complicated  by  the  ferromagnetic  nature 
of  the  steel.  Ferromagnetic  materials  are  characterized  by 

a.  Relatively  large  values  of  permeability 

b.  Variation  of  permeability  with  field  strength 

c.  Saturation  at  high  magnetic  field  intensities 


^’S.  S.  Brown  and  F.  W.  Anney,  Electrical  Resistivity  of  Line  Pipe  Steels, 
Paper  No.  69,  20th  Annual  Conference  of  National  Association  of  Cor- 
rosion Engineers  (1964). 
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Chemical  Composition  of  Line  Pipe  Steel  Specimens 


^ 00 
o o o o o o 
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*Five  samples  of  this  composition,  designated  by  6A,  6B,  6C,  60,  and  6E,  were  made  for  the 
electrical  resistivity  determinations.  (From  S.  S.  Brown  and  F.  W.  Anney,  Eleatvical 
Resistivity  of  Line  Pipe  Steels,  Paper  No.  69,  20th  Annual  Conference  of  National  Associa- 
tion of  Corrosion  Engineers  [1964].) 


Table  E5 


Comparison  of  Measured  and  Calculated  Resistivity* 
(Using  the  Percentage  of  Mn,  Si,  and  P and  the  Equation 
p = 11.59  + 5.43  [%Mn]  + 16.1  [%Si]  + 15.4  [%P]) 


Specimen 

No. 

• w,.#  ^ . -T  sy  L 'O'  " ’ J 

Measured 

Resistivity 

micro-ohm-cm 

' 1 V • 1 L 'O'-'  ' J 

Calculated 

Resistivity 

micro-ohm-cm 

* J / 

Deviation 

Measured 

Minus 

Calculated 

1 

15.51 

15.02 

.49 

2 

14.97 

14.98 

-.01 

3 

17.00 

17.04 

-.04 

4 

19.69 

19.66 

.03 

5 

20.05 

19.46 

.59 

6 

16.87 

17.54 

-.67 

7 

15.98 

16.73 

-.75 

8 

16.38 

16.63 

-.25 

9 

17.19 

17.10 

.09 

10 

19.77 

19.10 

.67 

11 

19.35 

19.62 

-.27 

12 

19.13 

1 9.36 

-.23 

13 

16.80 

16.71 

.09 

14 

17.24 

16.81 

.43 

15 

16.70 

16.91 

-.21 

16 

19.30 

19.56 

-.26 

17 

19.33 

19.03 

.30 

18 

19.56 

19.43 

.13 

19 

22.83 

22.80 

.03 

20 

21.39 

21.02 

.37 

21 

21.02 

21.51 

-.49 

22 

15.79 

1 5.58 

.21 

23 

15.66 

15.79 

-.13 

24 

15.26 

15.27 

-.01 

*From  S.  S.  Brown  and  F.  W.  Anney,  Eleatvioal  Resistivity  of  Line  Pipe 
Steels,  Paper  No.  69,  20th  Annual  Conference  of  National  Association 
of  Corrosion  Engineers  (1964). 
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d.  Hysteresis  (a  nonlinear  dependence  of  the  magnetization  on 
the  previous  magnetic  history  of  the  material). 

The  magnetic  properties  of  steel  are  often  strongly  dependent  on 
chemical  composition,  mechanical  working  and  heat  treatment.  Table 
E6  shows  nominal  values  of  the  relative  permeability 


for  several  iron  alloys  where 

pQ  = permeability  of  freo  space 
M = magnetization 
H = magnetic  field  intensity. 

It  is  apparent  that  magnetic  properties  are  greatly  influenced  by 
chemical  composition. 

Based  on  DC  resistance  measurements  and  RFI  attenuation  measure- 
ments, the  conductivity  and  the  single  value  of  relative  permeability 
u shown  in  Table  E7  have  been  suggested  for  the  purpose  of  calculating 
theoretical  shielding  effectiveness.  The  product  ya  based  on  these 
suggested  values  is  also  shown  in  Table  E7. 

Unfortunately,  the  permeability  of  conduit  steel  does  not  appear 
to  be  tabulated  in  the  literature. 


H 


Table  E6 


Relative  Permeability  of  Various  Ferromagnetic  Materials 


Material 

Initial 

Maximum 

Reference 

Very  pure  iron 

4,000 

8,000 

★ 

Silicon  iron 

3,500 

7,000 

Transformer  iron 

3,000 

5,500 

Machine  steel 

300 

450 

Cast  iron 

60 

90 

Low  carbon  iron 

200 

2,200-5,500 

1010  steel 

200 

3,800 

1%  Si  silicon  steel 

400-650 

1 ,700-6,000 

*Referenae  Data  fop  Radio  Engineers , 6th  Ed.  (H.  W.  Sams  & Co.,  1975), 
p 4-32. 

**T.  Lyman,  Editor,  Metals  Handbook,  8th  Ed.,  Vol  1 (American  Society  for 
Metals,  1966). 
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Table  E7 

Properties  of  Shielding  Materials  Based  on  DC  Resistance 
Measurements  and  RFI  Attenuation  Measurements* 


Material  Type 

Conductivity,  o 
(mho/m) 

Rel ati ve 
Permeabi 1 ity 

iLO. 

Galvanized  steel 

1.056X10^ 

270 

3580 

Hot  rolled  steel 

0.958X10^ 

160 

1930 

Cold  rolled  steel 

0.942X10^ 

127 

1500 

Terne  steel 

0.907X10^ 

157 

1790 

*From  R.  B.  Cowdell,  R.  A.  Hupp,  and  J.  N.  O'Leary,  RFI  Attenuating 
Materials  and  Structures , Technical  Report  AFAPL-TR-69-89  (Air  Force 
Aero  Propulsion  Laboratory,  Air  Force  Systems  Command,  1969),  p 125. 
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APPENDIX  F: 


EMP  PENETRATION  OF  DEFECTS,  CONDUIT 
FITTINGS,  AND  RELATED  HARDWARE 


Approach 

As  previously  indicated,  the  solid  conduit  provides  the  basic 
shielding  in  a conduit  system.  Defects  arising  in  the  conduit  such  as 
cracks  and  breaks  may  seriously  compromise  the  shielding.  Certain 
hardware  items  may  also  compromise  the  EMP  shielding  of  the  system. 
Completed  conduit  systems  may  include  numerous  fittings  and  related 
hardware  such  as  couplings,  unions,  case  access  fittings  (condulets), 
pull  boxes,  junction  boxes,  and  flexible  sections.  The  fittings  and 
hardware  used  are  usually  standard  conduit  accessories,  since  very  few 
fittings  and  hardware  items  have  been  designed  specifically  for  EMP- 
hardened  structures.  Thus,  these  items  may  compromise  the  EMP  shield- 
ing. In  order  to  assess  the  EMP  penetration  of  the  conduit  system,  it 
is  therefore  necessary  to  evaluate  the  leakage  contribution  as  a re- 
sult of  various  defects  and  hardware  items  for  the  current  waveshape  of 
interest. 

The  problem  of  determining  the  electromagnetic  field  leakage 
through  conduit  fittings  or  defects  is  less  amenable  to  direct  analy- 
tical solutions  than  is  the  problem  of  determining  the  electromagnetic 
field  penetration  of  solid  conduit.  The  analysis  is  complicated  by 
the  complex,  varied,  and  often  irregular  geometries  associated  with 
fittings  and  defects. 

Since  the  exact  mathematical  analysis  of  even  relatively  simple 
defects,  conduit  fittings,  and  hardware  items  is  a very  difficult  (and 
in  some  cases,  nearly  impossible)  task,  some  reliance  must  necessarily 
be  placed  on  empirical  results.  On  the  other  hand,  the  signals  re- 
sulting from  defects  quite  often  depend  on  the  magnitude  and  wave- 
shape of  the  applied  current  as  well  as  the  characteristics  of  the  de- 
fect itself.  Thus,  complete  reliance  on  empirical  data  would  necessi- 
tate acquisition  of  enormous  amounts  of  data  to  account  for  the  poten- 
tial variations  in  waveshape  and  characteristics  of  defects  which 
might  occur  in  actual  conduit  systems.  The  resultant  volume  of  data 
would  be  difficult  (if  not  impossible)  to  use,  and  extrapolation  to 
situations  for  which  data  were  not  taken  would  not  be  easy. 

Obviously,  some  compromise  between  a purely  theoretical  approach 
and  a purely  empirical  approach  is  required  to  arrive  at  a useful  re- 
sult for  the  designer  of  conduit  systems.  This  study  therefore  based 
its  approach  on  relatively  simple  mathematical  models  supplemented  with 
experimental  evaluation  of  certain  parameters. 
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The  signals  induced  on  an  internal  wire  within  a defective  con- 
duit can  be  considered  to  arise  as  a result  of  the  wire  coupling  with 
the  electric  and  magnetic  fields  associated  with  the  current  in  the 
vicinity  of  the  defect.  In  some  cases,  the  induced  signal  might  be 
interpreted  as  a result  of  direct  electric  field  and  magnetic  flux 
leakage.  Alternatively,  the  signal  might  be  considered  to  result  from 
the  impedance  of  a defect  changing  as  the  current  distribution  changes 
during  a pulse.  Or,  the  signals  might  be  interpreted  as  arising  from 
the  change  in  electric  and  magnetic  fields  at  the  outer  surface  of  a 
conduit,  since  the  electromagnetic  fields  would  be  expected  to  change 
as  the  current--which  initially  flows  almost  exclusively  on  the 
surface--penetrates  into  the  conduit  material.  Regardless  of  the 
exact  nature  of  their  origin,  the  observed  phenomena  can  be  described 
to  a first  approximation  by  a general  model  in  which  the  driving  vol- 
tage induced  by  a defect  has  a waveshape  that  is  proportional  to  the 
applied  conduit  current  plus  its  time  derivative.  This  appendix  con- 
siders some  details  of  this  general  model  for  EMP  penetrations  of 
defects. 


T-'. 
'<  11 


Electrical  Discontinuities 

Time  Domain 

To  a first  approximation,  the  electric  field  in  the  vicinity  of  a 
defect  may  be  considered  to  be  related  to  the  potential  across  a resis- 
tive component  of  the  defect  while  the  current  is  being  applied.  In 
other  words,  the  voltage  on  the  sense  wire  might  be  modeled  as 


VL(t)  = 


R^Ut) 


[Eg  Fi] 


where  R 


a coefficient  representing  the  linkage  between  the  conduit 
current  and  the  sense  wire. 


In  the  proposed  model  for  detects  for  a particular  position  of  the 
sense  wire,  the  coefficient  R-)  (in  units  of  ohms)  would  be  characteris- 
tic of  a given  defect  and  should  be  constant  for  that  defect.  Qualita- 
tively, one  would  expect  this  coefficient  to  vary  from  derect  to  defect 
depending  on  resistance  across  the  defect  and  proximity  of  the  sense 
wire  to  the  defect. 


have 


A contribution  to  the  leakage  signal  due  to 
the  same  waveshape  as  the  applied  current. 


this  effect  should 
In  other  words. 


‘o<‘) 

T^ty 


a constant 
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[Eg  F2] 


Frequency  Domain 

In  the  frequency  domain,  the  Fourier  transform  of  the  induced 
voltage 


V (t)  = I{t) 


[Eq  F3] 


V (oj)  = Koj) 


[Eq  F4] 


The  flaw  impedance  for  resistive  flaws  will  be  essentially  in- 
dependent of  frequency.  Thus,  the  conduit  current  and  the  conductor 
current  will  have  identical  waveforms: 


~ ® constant 


[Eq  F5] 


Rusty  couplings  are  an  example  of  electrical  discontinuities  which 
are  almost  entirely  resistive. 


Apertures 

Theoretical  Analysis 

Time  Domain.  In  the  presence  of  an  aperture,  a second  source  of 
a leakage  signal  could  be  linkage  of  magnetic  flux  with  the  sense  wire. 
Qualitatively,  one  would  expect  the  induced  voltage  on  the  sense  wire 
to  be  proportional  to  the  time  rate  of  change  of  the  flux,  i.e., 


V oc  . ^ 

V dt 


The  flux  is  proportional  to  the  magnetic  flux  density  B 


$ a:  B 


In  nonmagnetic  materials 


[Eq  F6] 


[Eq  F7] 


I 


It  follows  that  a contribution  from  this  source  would  be  expected  to 
be  of  the  form 


[Eq  FIO] 


where  = a coefficient  representing  the  linkage  between  the  magnetic 
flux  and  the  sense  wire 

In  the  proposed  model  for  defects,  the  coefficient  M2  (in  units 
of  henrys)  would  be  constant  for  a given  sense  wire  and  defect  confi- 
guration. However,  one  would  expect  the  coefficient  to  vary  with  de- 
fect dimensions  as  well  as  the  proximity  of  the  sense  wire  to  the 
defect.  A contribution  to  the  leakage  signal  from  this  mechanism 
would  be  expected  to  be  proportional  to  the  time  derivative  of  the  ap- 
plied current. 

The  general  model  for  the  direct  field  penetration  of  a defect  is 
therefore  proposed  to  be: 


V 


L 


R.|I  + M 


2 dt 


[Eq  Fll] 


Thus, 


VlU) 


[Eq  F14] 


Hence,  one  would  expect  the  transfer  impedance  to  be  frequency- 
dependent. 

HDL-measured  values  for  flaw  impedance  of  apertures  show  an  ap- 
proximate 20  dB  per  decade  increase  in  Zp.  This  implies  that 


2p  — cjci) 


[Eq  F15] 


where  c = a constant 
j = 

0)  = 2lTf 
f = frequency 

Thus,  from  Figure  FI  and  since  ^^<<2^  in  most  cases 


^0  = 


^0  ^0 


[Eq  F16] 


■ 2^  2w  du)  = 2^^  [^^ 

0 •'  -00  0 ^ 


The  constant 


c I 


can  be  evaluated  from  the  measured  flaw  impedance  curve. 


[Eq  F18] 
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mpedance  concept. 


It  should  be  noted  that  in  general  the  flaw  impedance  will  not  be 
an  easily  recognizable  function,  and  calculation  of  the  Fourier  trans- 
form of  the  function  will  be  more  complicated  than  this  example. 

Empirical  Results 

Time  Domain.  As  stated  previously,  the  matched  sense  wire  signal 
for  a conduit  having  an  aperture  has  been  hypothesized  to  include 
three  basic  components:  (a)  a small  contribution  due  to  diffusion,  (b) 
a resistive  component  proportional  to  the  current  applied  to  the  con- 
duit, and  (c)  the  component  due  to  induced  signals  resulting  from 
leakage  of  the  electromagnetic  fields  through  the  aperture.  The  third 
--or  leakage— component  is  further  hypothesized  to  be  proportional  to 
the  time  derivative  of  the  applied  conduit  current.  To  test  the  hypo- 
thetical model  for  the  matched  sense  wire  for  a conduit  with  an 
aperture,  a series  of  measurements  was  made  on  aluminum  conduit  samples. 
A solid  conduit  was  tested  to  show  the  diffusion  and  resistive  compon- 
ents of  the  signal.  Six  other  samples  containing  accurately  machined 
slots  1/2,  1,  2,  4,  8,  and  16  cm  long  were  tested  to  show  the  leakage 
signal  contribution.  The  slots  were  machined  to  the  lengths  stated 
with  a depth  of  exactly  one-half  of  the  conduit  diameter.  Thus,  for 
the  slot  length,  one-half  of  the  wall  was  machined  out.  Figure  F2 
shows  the  conduit  samples. 

The  sense  wire  signal  was  measured  with  the  wire  held  exactly 
centered  in  the  conduit.  Centering  was  accomplished  by  means  of  two 
plastic  (Lucite)  washers  placed  inside  the  conduit  with  a center  hole 
in  each  washer  for  passage  of  the  sense  wire. 

The  current  source  for  these  experiments  was  the  battery  dis- 
charge configuration  described  in  Appendix  B.  The  current  was  injected 
into  the  conduit  samples  by  means  of  12  AWG  copper  wires  attached  to 
the  conduit  using  stainless  steel  radiator  hose  clamps  (Figure  F3). 

To  minimize  the  resistive  and  diffusion  components  of  the  signal,  the 
current  was  injected  1 cm  from  the  slot  ends.  The  injection  wires  were 
placed  on  the  side  away  from  the  <;lot  ooenino  to  minimize  field  dis- 
tortion due  to  the  wires. 

Verification  of  the  mathematical  model  was  accomplished  using 
the  following  steps: 

a.  Measu’^ement  and  recording  of  the  injected  current  and  matched 
sense  wire  signals  for  the  six  conduit  slot  samples 

b.  Determination  of  the  constants  and  coefficients  of  the  equa- 
tion from  the  ERDAC  III  data 

c.  Plotting  of  the  ERDAC  III  data  using  a PDP-11  minicomputer, 
tape  data  transfer  system,  and  a Calcomp  565  plotter 
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d.  Superimposing  plots  of  the  mathematical  model  expression 
(with  constants  as  determined  in  step  b)  on  the  ERDAC  III  data  plots. 

Step  a.  A conduit  sample  having  an  8-cm-long  slot  with  injection 
1 cm  from  each  slot  end  (or  10-cm  distance  between  injection  points) 
was  compared  with  a solid  conduit  with  current  injection  over  a 10-cm 
distance.  Figure  F4  shows  the  ERDAC  III  signals  for  the  solid  conduit 
and  Figure  F5  shows  the  ERDAC  III  signals  for  the  conduit  sample  with 
an  8-cm  slot.  For  both  conduit  samples,  the  measured  current  peak  was 
369  ERDAC  III  units,  which  converts  to  227.4  A.  The  measured  sense 
wire  signal  for  the  solid  conduit  had  a peak  value  of  19  ERDAC  III 
units,  which  converts  to  1.86  mV.  For  the  8-cm  slot,  the  sense  wire 
signal  tail  was  asymptotic  to  a value  of  42.4  ERDAC  III  units,  which 
converts  to  4.14  mV. 

From  purely  resistive  considerations,  the  expected  sense  wire 
voltage  would  be  defined  by 


'^R  ^^c  ^ 2 


[Eq  FI  9] 


where  Vp  = the  voltage  across  the  sense  wire  terminating  resistor 

I^  = the  injected  conduit  current 

R|^  = the  resistance  per  unit  of  length 

L = the  conduit  length  over  which  current  is  injected 

The  factor  of  1/2  is  used  because  of  the  voltage  division  by  the  two 
sense  wire  matching  resistors. 

Through  measurement  of  resistivity,  the  aluminum  conduit  was 
found  to  have  a resistance  of  2.4  x 10'^  ohms/m.  Using  this  value  and 
the  227. 4-A  conduit  current,  the  anticipated  voltage  across  the  sense 
wire  terminating  resistor  for  a 10-cm  solid  conduit  is  2.73  mV.  With 
a sample  of  the  same  length  containing  an  8-cm  slot  the  same  wire 
voltage  should  be  4.91  mV.  The  difference  between  the  expected  and 
measured  values  may  be  partially  attributable  to  the  accuracy  of  the 
ERDAC  III  (0.5  percent  of  full  scale,  and  partially  to  the  method  of  in- 
jecting current,  wherein  current  density  within  the  10-cm  length  is  not 
uniform  within  the  conduit.  At  any  rate,  the  values  presented  do  show 
that  the  contribution  due  to  diffusion  current  in  the  time  frame  of 
interest  and  for  the  test  method  used  is  small. 

Step  b.  In  the  present  case,  the  applied  current  has  the  form 
shown  in  Eq  F20  and  Figure  F6: 


Channel  1 (upper)  0.1  V full  scale 
Channel  2 (lower)  1.0  V full  scale 
(Time  displayed:  720  psec) 

Figure  F4.  Solid  conduit,  ERDAC  III  recording  of  injected  current 
(lower  trace)  and  sense  wire  signal  (upper  trace). 


Channel  1 (upper)  0.1  V full  scale 
Channel  2 (lower)  1.0  V full  scale 
(Time  displayed:  720  psec) 


Figure  F5.  Conduit  with  8-cm  slot,  ERDAC  III  recording  of  injected 
current  (lower  trace)  and  sense  wire  signal  (upper  trace) 


I(t)  = Iq 


[Eq  E20] 


Figure  F6.  Form  of  applied  current. 

The  derivative  of  the  applied  current  has  the  form  shown  in 
Eq  F21  and  Figure  F7. 


di 

dt 


[Eq  F21] 


Figure  F7.  Form  of  the  derivative  of  the  applied  current. 

According  to  this  general  model,  the  total  leakage  signal  measured  will 
be 

^ ^0 
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Thus,  the  leakage  signal  will  have  the  typical  form  shown  in  Figure  F8. 


i 


[Eq  F26] 


Vl(o) 

W 

Thus,  if  the  measured  current  can  be  fit  to  the  form 

I(t)  = [l-e'^^^°]  [Eq  F27] 


by  appropriate  selection  of  the  parameters  I and  Tq.  then  the  coef- 
ficients R-j  and  can  be  evaluated  with  relative  ease. 

The  parameters  and  Tq  can  be  readily  deduced  from  the  applied 
current  by  noting  that 


and  that  when 


the  time  is 


I(t)  = I 

t-«o  0 


[Eq  F28] 


= l-e"''  = 0.632  [Eq  F293 


[Eq  F30] 


In  practice,  however,  there  is  some  variability  in  the  digitized 
data  for  a variety  of  reasons.  For  example,  when  monitoring  both  the 
applied  current  and  the  leakage  signal,  the  transient  digitizer  alter- 
nately samples  the  signals  at  a sampling  interval  of  1 psec.  Hence, 
data  points  for  a specific  signal  are  spaced  at  intervals  of  2 ysec, 
which  is  significant  compared  to  the  time  constants  of  the  excitation 
current  and  leakage  signals.  This  causes  some  difficulty  in  determin- 
ing exactly  the  time  of  the  start  of  the  current  pulse  as  well  as  the 
peak  value  of  the  leakage  signal.  This  situation  requires  a slight 
modification  in  the  fundamental  technique  of  evaluating  the  coefficients 
and  M2. 

Consider  the  applied  current 


] 


[Eq  F27] 
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F 


for  early  times, 
yields 


Note  that  for  small  values  of  x,  a series  expansion 


2 3 

„-x  _ T „ 4.  X . X . 

e -1-x  + jr+p"'*’ 


[Eq  F31] 


Hence 


and 


I(t)  = [1  - (1  - t/T^  + ....)] 


I. 


[Eq  F32] 


Kt)  = t 


Thus,  for  early  times,  the  current  pulse  is  very  nearly  linear. 


[Eq  F33] 


Similarly,  for  small  times,  series  expansion  of  the  exponentials 
yields 

VL(t)  ='^iIo[l  - (l-t/  0*-  • ♦ -*°[1  - t/TQ+.  . .]  [Eq  F34] 


so  that 


1 “1  HI 

V.  (t)  ' “ t • - ? ? 

0 0 0 


[Eq  F35] 


or 


In  ^7 

VL(t)  = ^ ^ [R^  > t 

0 0 0 


[Eq  F36] 


Hence,  for  very  early  times,  the  leakage  signal  should  also  be  very 
nearly  linear. 

The  first  few  data  points  for  the  applied  current  are  fit  to  a 
straight  line; 
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I(t)  = + B^t 


[Eq  F37] 


Similarly,  the  first  few  data  points  for  the  leakage  signal  are  fit 
to  a straight  line: 


VL(t)  = + B2t 


[Eq  F38] 


Then  note  that 


so  that 


A,  - H — 
2 2 


[Eq  F39] 


[Eq  F40] 


[Eq  F41] 


The  first  few  data  points  of  the  ERDAC  III  were  fit  to  a straight 
line  defined  by 


I(t)  = A^  + B^t 


[Eq  F37] 


This  straight-line  fit  was  accomplished  by  performing  a least  squares 
regression  to  a straight  line  using  a programmable  desk  top  calculator 
(HP-65)  which  yielded  the  coefficients  A]  and  B-]  in  terms  of  ERDAC  III 
units.  Using  these  coefficients  and  the  value  of  the  baseline  of  the 
conduit  current  I(t),  a value  of  the  start  time  T^  was  determined  from 


T = ^baseline  ~^1 


[Eq  F42] 


This  value  for  Tq  is  more  accurate  than  that  which  could  have  been  ob- 
tained by  using  the  first  ERDAC  III  point  above  baseline. 
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Similarly,  the  first  few  data  points  of  the  measured  leakage  sig- 
nal were  fit  to  a straight  line  defined  by 


VL(t)  = A2  + B2t  [Eq  F38] 


with  A2  and  62  calculated  using  the  HP-65.  Using  T (obtained  above) 
in  this  equation,  a value  for  the  peak  leakage  signal  (Iq/to)  was  then 
determined  from 


''lI")  = ' ''2  " Vc 


[EQ  F43] 


The  value  for  was  determined  from  the  ERDAC  III  data  by  aver- 

aging the  baseline  of  the  injected-current  waveform  and  subtracting 
this  value  from  the  averaged  value  to  which  the  injected  current  rises 
asymptotically.  The  value  thus  obtained  was  then  converted  back  to 
amperes  from  the  ERDAC  III  units.  The  value  for  Tq  was  then  determined 
by  interpolation  between  ERDAC  III  data  points  to  determine  the  time 
when  0.632  (1/e)  times  the  final  value  is  reached  relative  to  the  value 
of  Tq  previously  determined. 

Steps  a and  d.  A PDP-11  minicomputer,  tape  data  transfer  unit, 
and  Calcomp  565  plotter  were  used  to  assist  in  verifying  the  mathemat- 
ical model  for  the  leakage  signal  through  the  conduit  slots.  The  ERDAC 
III  records  which  were  recorded  on  tape  were  transferred  to  the  tape 
data  transfer  unit,  and  from  there  were  input  into  the  PDP-11  computer. 
The  PDP-11  was  programmed  to  plot  the  waveforms  as  recorded  by  the  ERDAC. 
In  addition,  a program  was  written  to  enable  the  PDP-11  and  Calcomp 
plotter  to  compute  and  plot  a summation  of  exponentials  of  the  form 

f(t)  = ae““^  + be“®^  + ce"^^  + de"*^^  + ge"®^  [Eq  F44] 


Thus,  after  plotting  the  conduit  current  waveform 

- t/T 

I(t)  = Iq  (l-e  0)  [Eq  F27] 

the  ERDAC  III  data  were  used  to  derive  values  for  Iq  (the  current  maximum 
value)  and  Tq  (the  1/2  or  e-fold  value).  In  f(t)  above  Uq  F44),  c,  d, 
and  g were  set  to  zero,  a and  b were  set  equal  to  Iq,  a was  set  to  in- 
finity, and  8 was  set  equal  to  1/tq.  This  curve  was  then  plotted  as  a 
continuous  line  superimposed  onto  the  plot  of  the  ERDAC  recorded  wave- 
form, which  was  represented  by  a sequence  of  squares.  Figures  F9,  FIO, 
Fll,  and  F12  show  these  curves  for  the  noted  slot  lengths.  Excellent 
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agreement  is  shown  between  the  calculated  curve  and  the  experimental 
data  curve. 

The  leakage  signal  has  been  hypothesized  to  be  of  the  form 


-t/x  I -t/x 

VL(t)  = RiIq  - Riloe  +^2:^0  ° [Eq  F45] 

^0 


Thus,  in  Eq  F44,  d and  g were  set  to  zero,  a was  set  to  infinity,  and 
6 and  y were  set  equal  to  Xq.  A value  for  R-j  Iq  was  derived  from  the 
ERDAC  III  data  as  the  asymptotic  value  of  V|^  and  the  value  for 


was  derived  from  V|_(o),  the  peak  value  of  Vl.  Plots  of  the  ERDAC  re- 
corded waveforms  were  made  by  a series  of  squares,  and  the  calculated 
waveforms  were  superimposed  as  continuous  lines.  Figures  F13,  F14,  F15 
and  F16  show  the  results  for  the  2,  4,  8,  and  16-cm  slots.  Plots  were 
not  made  for  the  1/2-  and  1-cm  slots  because  of  the  low  signal  magni- 
tudes. Each  plot  shows  reasonably  close  agreement  between  the  calcu- 
lated and  experimental  curves.  The  conclusion  can  thus  be  drawn,  at 
least  to  a first  approximation,  that  the  proposed  model  does  define  the 
sense  wire  signal. 

It  is  therefore  seen  that  the  mathematical  expression  and  a table 
of  values  for  the  coefficients  R-]  and  M2  can  be  used  to  define  every 
slot  type  of  defect  in  a conduit.  An  hypothesis  at  this  point  is  that 
a similar  model  will  apply  to  more  complex  apertures,  such  as  loose 
coupling  threads,  unions,  and  other  conduit  fittings. 

Table  FI  provides  values  for  R]  and  M2  for  the  six  slots  evaluated. 
Each  value  listed  represents  an  average  of  two  trials. 

Table  FI 

Values  of  R^  and  M2 


Slot  Size,  cm 

R^ , ohms 

M^,  henries 

1/2 

5.8x10"® 

3.04x10"^® 

1 

6.36x10"® 

5.90x10"^® 

2 

7.44x10"® 

8.96x10'^° 

4 

1.21x10'® 

17.1x10"^® 

8 

1 .80x10"® 

37.6x10"^® 

16 

2.99x10'® 

63.8x10"’® 
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An  additional  objective  of  the  empirical  rk  was  to  determine 
how  the  peak  leakage  signal  for  a given  conduit  current  rise  time  var- 
ies with  slot  length.  Thus,  the  problem  was  that  of  determining  the 
coefficient  M2  versus  slot  length.  In  order  to  make  this  determination, 
the  six  aluminum  conduit  samples  with  slot  lengths  of  1/2,  1,  2,  4,  8, 
and  16  cm  were  tested.  Two  sets  of  data  were  taken  for  each  sample, 
using  the  ERDAC  III.  For  each  data  set,  the  data  interpretation 
analysis  described  on  p 99  was  performed.  The  data  derived  are  sum- 
marized in  Table  F2  and  plotted  in  Figure  F17.  This  plot  shows  an 
approximate  straight-line  relationship  between  M2  and  slot  length  for 
the  larger  slot  lengths.  Extrapolation  of  the  data  points,  however, 
indicates  that  the  curve  would  not  pass  through  the  origin. 

Table  F2 

Relationship  of  M and  Slot  Length 


Sample 

Trial  1 

Vl(0)/B^ 

Trial  2 

Vl(0)/B^ 

Average 

Vl(0)/B^ 

2-,  1/2-cm  slot 

.91 

1.01 

.96x9. 77x10"^ 

3-,  1-cm  slot 

1.95 

1.78 

1.87x9.77x10"^ 

4-,  2-cm  slot 

3.01 

2.66 

2.84x9.77x10".^ 

• * 

5-,  4-cm  slot 

5.16 

5.64 

5.40x9.77x10"^ 

6-,  8-cm  slot 

11.49 

12.31 

11.90x9.77x10"^ 

7-,  16-cm  slot 

20.00 

20.30 

20.15x9.77x10"^ 

Flaw  Impedance  of  Apertures  (Transverse  Slots).  HDL  measured  the 
flaw  impedance  of  transverse  slots.  Figure  FI 8 shows  the  flaw  impedance 
of  a transverse  slot  in  a 4-in.  (102-mm)  conduit  with  the  internal  con- 
ductor (sense  wire)  centered  in  the  conduit.  The  location  of  the 
internal  conductor  in  relation  to  the  slot  has  a considerable  effect  on 
the  flaw  impedance,  as  shown  in  Figure  F19.  Figure  F20  shows  the  worst 
case  for  a slot  for  which  the  flaw  impedance  was  measured.  This  worst 
case  was  a flexible  conduit  such  as  those  used  in  the  SAFEGUARD  BMD 
system  completely  removed  leaving  only  a parallel  copper  strap  connect- 
ing the  conduit  sections. 
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RELATIVE  V,  (0)/B, 


WIDTH  = 1/2 


gure  F18.  Flaw  impedance  of  a transverse  slot  in  a 4-in.  (102-mm)  inner  diameter  conduit 
with  the  sense  wire  centered  in  the  conduit.  (From  H.  A.  Roberts,  J.  Capobianco 
and  F.  Agee,  SAFEGUARD  Buried  Conduit  Studies  [HDL,  undated]). 


Figure  FI  9.  Flaw  impedance  of  transverse  slot  with  sense  wire  at 
three  different  positions.  (From  H.  A.  Roberts,  J. 
Capobianco,  and  F.  Agee,  SAFEGUARD  Buried  Conduit 
Studies  [HDL,  undated]). 


FREQUENCY  (MHi) 

Figure  F20.  Flaw  impedance  of  large  aperture.  (From  H.  A.  Roberts, 
J.  Capobianco,  and  F.  Agee,  SAFEGUARD  Buried  Conduit 
Studies  [HDL,  undated]). 
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APPENDIX  G: 
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j COEFFICIENTS  FOR  CALCULATION  OF  EMP 

PENETRATION  OF  DEFECTS,  CONDUIT 
j FITTINGS^  AND  RELATED  HARDWARE 

I 

I This  appendix  supplies  coefficients  for  the  mathematical  model  to 

i enable  a simple  prediction  of  the  anticipated  EMP-induced  disturbance 

I for  a number  of  types  of  conduit-related  fittings  or  hardware  items. 

! Results  of  a time  domain  investigation  of  1-in.  (25-niii)  aluminum 

i couplings,  unions,  and  a type  C condulet  (case  access  fitting),  and  a 

j frequency  domain  flaw  impedance  investigation  of  rusted  couplings  and 

i unions  are  summarized. 


Couplings 

Time  Domain 


i > ' 


1 


An  aluminum  coupling  (clean  threads)  was  tested  to  determine  the 
leakage  and  resistive  components  on  a matched  sense  wire  for  varying 
degrees  of  tightness.  The  current  injection  method  used  was  identical 
to  the  method  used  in  testing  slots  (Appendix  F).  For  each  tightness 
level  tested,  the  resistance  of  the  coupling  and  thread  matings  was 
measured  by  the  DC  current  injection  and  voltage  drop  measurement 
method. 

Figure  G1  shows  the  recorded  ERDAC  III  waveform  for  injected  cur- 
rent and  the  matched  sense  wire  signal  for  a very  loose  coupling  (1/2 
turns  mated  on  each  end).  This  photograph  indicates  that  the  leakage 


Upper  trace  - sense  wire  voltage,  0.1  V full  scale 
Lower  trace  - conduit  current,  632  A full  scale 

Figure  G1 . ERDAC  III  waveforms  for  loose  coupling. 

no 


signal  is  too  small  to  measure  relative  to  the  resistive  component. 
Careful  examination  also  shows  that  the  rise  time  for  the  sense  wire 
signal  is  also  longer  than  the  rise  time  of  the  injected  current.  This 
delay  in  rise  time  is  attributable  to  the  diffusion  delay  time. 

Table  61  presents  peak  sense  wire  signal  versus  coupling  tightness 
and  resistance. 


Table  61 

Peak  Sense  Wire  Signals  vs.  Coupling 
Tightness  and  Resistance 


Tightness 

DC 

Resistance, 
ohms  X 10'^ 

Peak  Conduit 

Current  (I  ) 
c 

amps 

Peak  Sense 
Wire  Signals 
i\),  mV 

Hand  tight 
less  1 turn 

.307 

123.5 

109.0 

Hand  tight 

0.011 

121.0 

19.5 

Wrench  tight 

0.008 

120.4 

18.6 

It  should  be  noted  that  the  data  in  Table  61  were  derived  with 
current  injection  over  about  18  in.  (0.46  m)  of  conduit,  including  the 
coupling.  The  signal  measured  thus  includes  a substantial  contribution 
due  to  diffusion  and  resistive  signal  components. 

Conduit  Coupling  Flcoj  Impedance 

Figure  62  shows  the  HDL-measured  flaw  impedance  for  a rusted  coupling. 
The  shape  of  the  flaw  impedance  curve  will  probably  be  similar  for 
different  thicknesses  of  the  rust  layer;  however,  its  magnitude  will  be 
different.  The  impedance  is  virtually  independent  of  frequency  for  the 
rusty  coupling  joint  at  frequencies  less  than  10  MHz. 


Unions 

Time  Domain 

An  aluminum  explosion-proof  union  was  tested  as  described  in  the 
previous  section.  The  union  ends  were  wrench-tightened  onto  the  mating 
conduit  sections,  but  the  coupling  ring  was  set  at  varying  tightnesses 
for  the  tests.  In  all  tests  performed,  the  leakage  signal  was  too  small 
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to  measure.  The  resistive  component  and  diffusion  components  were  thus  | 

the  only  signal  portions  measured.  Since  injection  was  made  over  about 
18  in.  (0.46  m)  of  the  conduit  including  the  union,  the  constant  dif- 
fusion and  resistive  components  for  that  portion  of  the  sample  prevent  a 
correlation  of  sense  wire  signal  versus  the  measured  resistance  of  the 
union.  Table  G2  presents  data  from  the  tests. 

Conduit  Union  Flaw  Impedance 

HDL  measured  the  flaw  impedance  for  a 4-in.  (102-mm)  conduit  union 
for  five  cases:  (a)  a relatively  loose  assembly  at  50  ft-lb  (68  N-m) 
of  torque,  (b)  a tight  assembly  at  300  ft-lb  (407  N-m)  of  torque,  (c) 
an  assembly  having  a teflon  washer  placed  between  the  mating  faces  of 
the  union  pieces,  thus  preventing  contact,  (d)  an  assembly  having  teflon 
insulation  placed  between  concentric  parts  of  the  union,  preventing 
adequate  electrical  contact,  and  (e)  an  assembly  having  a Lucite  wedge 
placed  between  the  mating  faces  of  the  union,  preventing  even  contact. 

Figure  G3  presents  the  results  of  these  measurements. 


Table  G2 

Union  Time  Domain  Test  Data 


Tightness 

Conduit  Current, 

(Ij.)  amps 

Sense  Wire 

Signal , mV 

Resistance, 

mn 

Hand  tight 

2 turns 

227.3 

37.9 

0.26 

Hand  tight 

1/2  turn 

225.4 

10.6 

0.071 

Hand  tight 

1/8  turn 

231 .6 

6.2 

0.036 

Hand  tight 

231.6 

3.7 

0.022 

Wrench  tight 

242.1 

3.2 

0.018 

Conduit  Fittings 

An  aluminum  conduit  fitting  (Type  C,  straight-through,  1-in. 
[25.4-mm]  Electrolet*)  was  tested  as  described  for  the  couplings.  The 


*Electrolet  is  a brand  name  of  the  Killark  Corporation. 
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fitting  was  wrench-tightened  onto  the  mating  conduit  sections  at  each 
end.  Tests  were  made  with  and  without  the  aluminum  cover  plate  on. 
Figure  G4  shows  the  ERDAC  III  recorded  signals  for  the  fitting  with 
the  cover  off.  A leakage  component  exists,  but  it  is  much  smaller 
than  the  leakage  component  for  a conduit  with  slot  length  and  size 
comparable  to  the  opening  in  the  fitting.  This  difference  is  attribu- 
table to  the  fact  that  the  fitting  flange  extends  well  beyond  the 
sense  wire  and  thus  the  fields  must  penetrate  much  more  deeply  than  in 
the  slotted  sample. 

When  tested  with  the  cover  on,  the  leakage  signal  component  on 
the  sense  wire  became  too  small  to  measure,  but  the  diffusion  and 
resistive  signal  components  remained  essentially  the  same.  Table  G3 
presents  data  for  the  tests. 


Table  G3 

Conduit  Fitting  Time  Domain  Test  Data 


Peak  Conduit 

Peak  Leakage 

Peak  Resistive 

Current,  A 

Signal , mV 

Signal , mV 

Lid  off 

221.7 

3.3 

3.8 

Lid  on 

227.3 

0 

4.1 

i 

, 

^ 1 

Top  trace  = sense  wire  signal,  0.1  V full  scale 
Bottom  trace  = conduit  current,  632.4  A full  scale 
(time  base;  1 psec/sample,  725  psec  displayed) 

Figure  G4.  ERDAC  III  waveforms  for  Electrolet,  Type  C with  cover  off. 
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APPENDIX  H: 

EFFECTS  OF  CIRCUIT  CONFIGURATION  ON 
EMP- INDUCED  SIGNALS 


Transmission  Line  Propagation  of  Signals 

For  rapidly  varying  EMP-induced  signals,  the  propagation  must  be 
analyzed  using  transmission  line  theory.  The  onset  of  the  signal  at 
the  end  of  a conduit  system  will  be  delayed  by  the  time  it  takes  the 
signal  to  travel  the  distance  from  the  location  of  the  defect  to  the 
end  of  the  conduit.  Furthermore,  possible  reflections  due  to  mismatch 
in  termination  of  a transmission  line  may  also  arise. 


Lumped  Parameter  Circuits 

For  slowly  varying  signals,  the  transmission  line  properties  be- 
come less  significant,  and  in  the  limit  such  problems  approach  the 
equivalent  circuit  problem  with  lumped  parameters.  The  time  delay  due 
to  propagation  time  becomes  negligible  for  slowly  varying  signals, 
since  the  velocity  of  propagation  is  nearly  that  of  the  speed  of  light 
(c  = 3 X 10°  m/sec). 


EMP-Induced  Voltages  and  Currents 

Open  Circuit  Voltage 

Time  Domain  Theoretical  Analysis.  The  open  circuit  voltage  in- 
duced on  a sense  wire  within  a conduit  can  be  expressed  as  the  EMP- 
induced  electric  field  integrated  along  the  length  of  the  conduit. 

Eq  HI  gives  the  relation  for  this  voltage: 


E(a,e)  = JIQ[ 


yo2(b-a)^ 


-l2n4Tif 

[•re  ^ |-(2n+l)  1-11 


n=0  e 


[Eq  HI] 


where  % = the  length 

Q = the  injected  charge 
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a 


inner  radius 


b = outer  radius  ] 

I 

0 = conductivity  | 

u = permeability  I 


0 = a relative  time  expressed  in  units  of  a characteristic 
time  0 = t/T 


where 


■[■  _ ua 


Eq  HI  is  the  product  of  the  electric  field  expression  (Eq  C15)  derived 
in  Appendix  C and  the  length  of  the  conduit. 

Eq  HI  converges  well  for  short  times,  which  is  convenient  for 
generation  of  the  response  to  an  impulse  current  injected  on  the  surface 
of  the  conduit. 


Frequency  Domain.  To  derive  an  expression  for  the  open  circuit  vol- 
tage (Eq  HI),  the  EMP-induced  internal  electric  field  was  integrated 
along  the  length  of  the  conduit.  After  this  integration,  all  that  was 
introduced  to  modify  the  expression  for  the  internal  electric  field  was 
the  length  of  the  conduit.  Because  of  this  time- independent  modifica- 
tion, the  Fourier  transform  of  the  open-circuit  voltage  expression  is 
equivalent  to  the  Fourier  transform  of  the  internal  electric  field; 
therefore,  the  open  circuit  voltage  frequency  response  is  equivalent  to 
the  induced  electric  field  frequency  response,  which  is  investigated  in 
Appendix  C. 


Experimental  Evaluation.  To  determine  if  Eq  HI  is  accurate  in 
describing  the  amplitude  and  time  response  of  the  open-circuit  voltage 
induced  by  an  injected-current  pulse  which  approximates  an  impulse 


function,  the  experimental  procedure  described  in  the  Experimental 
Evaluation  section  of  Appendix  C was  used. 


A 3-m  length  of  nominal  1-in  (25-mm)  aluminum  conduit  was  placed 
in  the  setup  shown  in  Figure  C4.  A current  pulse  with  an  amplitude  of 
350  A,  a rise  time  of  10  nsec,  and  a fall  time  of  1.5  psec  was  injected 
onto  the  surface  of  the  conduit.  The  rise  and  fall  times  of  the  in- 
jected current  pulse  are  much  shorter  in  duration  than  any  of  the  time 
scales  of  the  response.  The  open-circuit  voltage  was  measured  with  a 
digital  recording  oscilloscope  which  gave  a high  degree  of  accuracy  for 
the  data  points. 


To  determine  the  experimental  agreement  of  Eq  HI  and  the  data 
points,  a plot  of  Eq  HI  was  fit  to  the  experimental  data  for  the  open 


circuit  voltage  by  selecting  values  for  T and  the  quantity 

r . ^ 1 

The  experimental  data  points  were  normalized  and  plotted  on  the  theore- 
tical curve  calculated  using  the  selected  values  for  the  parameters. 
Empirically,  it  was  found  that  the  theoretical  curve  was  in  excellent 
agreement  with  the  experimental  data  points  if  the  following  values  for 
the  parameters  were  selected: 


T = 85  psec 


T^^^/ab  pa^(b-a)^ 


} = 3.38  X lO'"’  volts 


As  can  be  seen  from  its  plot  (Figure  HI),  the  time  response  of 
Eq  HI  is  in  excellent  agreement  with  the  experimental  data. 

To  evaluate  the  accuracy  of  the  theoretical  relation  used  to  de- 
termine the  amplitude  of  the  open  circuit  voltage,  the  constant  in 
Eq  HI , 


i^Q{ 


pa^(b-a)^ 


was  calculated  using  parameters  of  the  conduit  and  the  injected  charge. 
This  value  was  then  compared  with  the  value  required  in  the  theoretical 
calculation  of  the  peak  amplitude  to  be  equal  to  the  experimental  peak 
ampl itude. 

As  previously  noted,  the  empirical  value  of  the  amplitude  parameter 
was  determined  to  be  3.38  x lO"^  V.  The  theoretical  value  was  calculated 
to  be  3.36  x 10"3  using: 

a = .0133  m 

b = 0.0167  m 

p = Att  X 10"^  henrys/m 

0 = 4 X 10"^  c 

a = 2.64  X 10^  mhos/m 

The  theoretical  value  is  within  1 percent  of  the  empirical  value. 
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In  summary,  the  theoretical  expression  for  the  open  circuit  vol- 
tage is  in  excellent  agreement  with  the  experimentally  determined  open 
circuit  voltage. 

Short  Circuit  Current — Theoretical  Analysis 

Equivalent  Circuit.  In  modeling  the  short  circuit  current  induced 
on  a sense  wire  inside  a conduit  which  has  an  injected  current  pulse,  a 
lumped  parameter  equivalent  circuit  can  describe  the  system.  Figure  H2 
shows  such  a circuit  consisting  of  a voltage  source  (which  is  the  open 
circuit  voltage  induced  on  the  sense  wire),  the  resistance  of  the  sense 
wire  and  conduit,  and  the  inductance  of  the  sense  wire-conduit  combina- 
tion. Capacitance  is  neglected  in  the  model. 

Writing  the  differential  equation  for  the  current  in  the  loop 

gives 


L ^ + RI  = (a,t) 


[Eq  H2] 


The  Laplace  transform  of  this  equation  is 


(Ls  + RT 


[Eq  H3] 


From  Eq  Cl  4, 


E (a,s)  = -9—  /S  [ 

^ 2iT»/ab  ^ sinh[(b-a)/uoa3 


-I 


[Eq  H4] 


Thereto  re. 


/s 


l(s)  = A]  • [—  ^ 

2TT/ab  sinh[(b-a)/uoaj 


Introducing  the  conduit  diffusion  time  T,  where  T = — 


^ ' ITTr  [Eq  H5] 
2 


i(s)  = ys 

/nr  / o 


2iT/ab  [sinh[2/r/s]][Ls  + R] 


[Eq  H6] 


1 


I 

I 
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where  I = the  length 


Q = the  charge  injected 
b = outer  radius 
a = inner  radius 
y = permeability 
a = conductivity 

To  return  to  the  time  domain,  the  inverse  Laplace  transform  of 
Eg  H6  was  taken.  The  method  used  here  to  invert  is  the  method  of 
residues,  which  is  described  in  the  annex  to  Appendix  C.  The  time 
domain  solution  for  the  short  circuit  current  is  given  by  Eq  H7 


I(t)  = A 


[Eq  H7] 


Eq  H7  is  the  short  circuit  current  response  to  an  impulse  current  pulse, 
in  terms  of  relative  time  6 (which  is  the  ratio  of  the  diffusion  time 
and  time)  and  a dimensionless  parameter  which  is  the  ratio  of  the  con- 
duit diffusion  time  T and  the  L/R  time  constant. 

Properties  of  Short  Circuit  Current  Response.  Starting  with  Eq  H7, 
it  is  possible  to  investigate  such  properties  as  where  in  time  the  peak 
current  occurs,  the  amplitude  of  the  peak  current,  the  relationship 
between  the  parameter  a and  the  peak  amplitude  and  the  time  of  the  peak, 
and  the  relationship  between  the  amplitude  of  the  peak  and  time  of  the 
peak.  These  properties  are  investigated  in  the  following  two  sections. 

Time  of  Peak  Current.  The  maximum  current  occurs  at  the  time  for 
which  ^ = 0.  That  is,  when 
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3/2  -a0  - (-i)"(.'yL_)n‘'  e 


TA  ? + I 

sin  (2/a)  n=l 


n 2 2, 

vO,  n IT  » 2 

) g 

P 2 4cx  -1 
[n  - -y] 


= 0 [Eq  H8,] 


Therefore,  the  time  of  the  peak  occurs  at  some  6 peak  such  that 


sin(2/a) 


/ oo 

peak  ^ ^ 
n=l 


4 2- 
(-1)"  -4—  e 

- —7] 


[Eq  H9] 


The  above  relation  involves  finding  the  zero  of  an  infinite  series, 
but  if  0pgg|^  is  ~1  or  larger,  a first  approximation  of  0pg3|^  can  be 

found  using  the  first  term  of  the  infinite  series.  The  first  approxima- 
tion could  then  be  used  to  determine  the  value  of  0_p,|.  by  iteration 
which  gives  a maximum  for  a specified  value  of  a.  ^ 

©peak  determined  on  a programmable  calculator  by  iteration. 

The  minimum  value  of  Qpg^j^  is  .36701  for  large  a.  This  is  also  the 

value  of  the  open  circuit  voltage  peak  (see  Appendix  C).  As  a 
approaches  zero,  the  time  at  which  the  peak  occurs  becomes  longer. 
Figures  H3  and  H4  are  plots  of  the  time  that  the  peak  occurs  (relative 
time  0)  versus  a. 

The  time  at  which  the  short  circuit  peak  will  occur  can  be  calcu- 
lated from  Figures  H3  and  H4  if  the  diffusion  time  of  the  conduit 
(Table  C2)  and  the  value  of  the  time  constant  L/R  are  known.  The  time 
of  the  peak  is  given  by  the  relation: 


^peak  ' ®peak  ^ 


[Eq.  HIO] 


Peak  Amjilitude  of  Cuvrent.  The  expression  in  the  brackets  in 
Eq  H7  was  evaluated  at  various  values  of  0pgg|^  determined  in  the  pre- 
vious section.  Figures  H5  and  H6  plot  the  peak  amplitude  of  -the  expres- 
sion in  the  brackets  versus  Qpggi^-  As  can  be  seen  from  the  plots,  the 

quantity  in  the  brackets  reaches  a limiting  value  of  .5  as  Qpg^i^  gets 

larger,  which  implies  a approaches  zero.  The  peak  amplitude  of  the  ex- 
pression in  the  brackets  was  also  plotted  versus  the  parameter  a 
(Figure  H7). 

The  amplitude  of  the  peak  current  can  be  determined  if  6pgg|^  is 

known  (Figures  H5  and  H6)  or  if  the  parameter  a is  known  (Figure  H7). 

The  peak  current  amplitude  is  then  given  by: 
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T(R/L) 


Figure  H6.  Plot  of  magnitude  peak  vs. 
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T = rAS.  1 _L  1 . peak  amplitude  of  term  [eq  nil] 
peak  '-2tt  / a L ^ ’ in  brackets 

I 

Short  Civauit  Current  Frequency — Domain 

Although  the  main  objective  of  this  study  is  to  examine  the  time 
variation  in  the  short  circuit  current,  it  is  useful  to  consider  the 
frequency  domain  of  the  short  circuit  current  transfer  function. 

Fourier  Transform.  From  the  definition  of  the  Fourier  transform, 
the  transfer  function  for  the  short  circuit  current  equation  (Eq  H6) 
can  be  modified  by  the  simple  substitution  of  s = ju  to  obtain  the 
Fourier  transform.  After  this  substitution,  Eq  H6  becomes 


' 2TT/ib  sinh(2/j^T)’ 

Multiplying  numerator  and  denominator  by  (R  - jwL)  and  letting  /T 
v^/2  (1  + j). 


[Eq  HI  2] 


jw) 


/^]  [-  j ^Eq  HI  3] 

2tt/^  ° R'^+io  L sinh(/^  + jv^) 


Using  transcendental  function  identities: 

sinh  (x  + y)  = sinh  x cosh  y + cosh  x sinh  y 
cosh  (jB)  = cos  B 
sinh  (jB)  = j sin  B 


H 


i 

1 


I = r_A-5_ 

2tt/^ R 


[- 


(R  - j(^L) 


/2^ 


s i n h ( ) c 0 s ( /2(irr)  + j s i n ( ) c 0 s h 


The  magnitude  of  the  transfer  function  is  given  by: 


[Eq  H14]  ; 

i 
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[Eq  HI  5] 


/ R^+aj^L^(sinh^/2u)T  cos^i/2ajT+sin^/5wT 


The  phase  of  the  transfer  function  is  given  by; 


({),.  V = tan  [ 


(l-ifjR)  - (1 +‘^)(coth/2;7r  tan^) 

" ' [ *- t ■ 

(1  + <i|R  + (1  . yB.)(cothAyr  tan  y^)' 


-]  [Eq  HI  6] 


When  the  relative  frequency  5 is  defined  as  3 wT  and  the  parameter  a 
as  TR/L,  the  magnitude  is  given  by 


|GU»)|  = /¥ 


27rv^iF 


[Eq  HI  7] 


2R^  cos^(^)  ^ sin^  (/p  cosh^  (/^jj 

o2  _ 2 \ ^ 2"  t / 


and  the  phase  by: 


= tan 


[1  -(|^)]-  [1  coth  tan  (V|) 

C'  cl”)  ta"  c?) 


[Eq  H18] 


Properties  of  Transfer  Function.  To  determine  the  properties  of 
the  transfer  function,  both  magnitude  and  phase  were  plotted.  Figure  H8 
is  a plot  of  the  magnitude  of  the  transfer  function  in  decibels  versus 
the  relative  frequency  C.  The  transfer  function  magnitude  was  plotted 
for  a range  of  0.1  to  °°  for  a.  As  can  be  seen  from  the  plot,  the 
magnitude  is  reduced  by  150  dB  within  a relative  frequency  range  of  800 
to  1100  for  the  range  of  a.  To  convert  the  relative  frequency  ^ to 
frequency,  the  relationship  f(Hz)  = ^/8T2it  can  be  employed.  For  1 - in. 
(25-mm)  aluminum  conduit,  T = 85  usee,  the  relative  frequency  range  of 
800  to  1100  converts  to  190  kHz  and  260  kHz.  Thus,  for  1-in.  (25-mm) 
aluminum  conduit  the  attenuation  is  150  dB  at  200  kHz,  and  at  1 MHz 
the  attenuation  is  greater  than  250  dB.  For  1-in.  (25-mm)  galvanized 
steel  conduit,  T = 4.2  msec  (Table  C2),  the  attenuation  is  150  dB  at 
about  4 kHz,  and  at  MHz  the  attenuation  is  greater  than  1200  dB. 


The  phase  of  the  transfer  function  was  plotted  for  a range  of  a 
of  0.1  to  and  a relative  frequency  range  of  1 to  1000  (Figure  H9). 

Short  Circuit  Current  ~ Experimental  Evaluation 

The  theoretical  expression  Eq  H7  was  verified  experimentally  by 
measuring  the  short  circuit  current  induced  on  the  sense  wire  with  a 
high  current  pulse  injected  on  the  surface  of  the  conduit.  The  conduit 
used  in  this  test  was  a 3-m  length  of  nominal  1-in.  (25-mm)  aluminum 
conduit. 

Figure  HIO  shows  the  experimental  setup.  A current  pulse  with  an 
amplitude  of  400  A,  a rise  time  of  10  nsec,  and  a fall  time  of  1 psec 
was  injected  on  the  surface  of  the  conduit.  The  current  pulse  was  such 
that  it  was  a good  approximation  of  an  impulse  function,  and  thus  pro- 
vided a means  of  directly  verifying  the  theoretical  expression. 

To  measure  the  current  induced  on  a wire  inside  the  conduit,  a 
sense  wire  was  inserted  into  the  test  conduit  as  shown  in  Figure  HIO. 
The  sense  wire  was  connected  to  the  pulser  end  of  the  conduit  through 
a resistor  with  the  end  of  the  conduit  capped  to  keep  leakage  signals 
from  entering  the  system.  The  sense  wire  ran  the  entire  length  of  the 
conduit  and  extended  into  the  shielded  room  where  it  was  shorted  to 
the  grounded  wall  of  the  shielded  room.  A Pearson  No.  411  inductive 
current  probe  was  placed  around  the  sense  wire.  Output  of  the  current 
probe  was  then  amplified  by  a high-speed  operational  amplifier  and 
measured  by  a digital  recording  oscilloscope.  Data  points  were  then 
accurately  read  from  the  trace  displayed  by  the  scope. 

Tests  were  also  conducted  with  additional  resistance  in  series 
with  the  sense  wire  inside  a shielded  room. 

The  resistance  of  the  entire  system  was  determined  by  injecting  a 
DC  current  and  measuring  the  voltage  drop.  The  inductance  of  the 
system  was  measured  using  a Tektronix  Model  No.  130  L-C  meter. 

To  determine  the  agreement  of  Eq  H7  to  the  data  points  measured, 
the  theoretical  curve  was  fit  to  the  experimental  data  for  the  short- 
circuit  current  by  selecting  values  for  T,  L/R,  and  the  quantity: 

i -1 

2-n/ab  /f 

The  experimental  data  points  were  normalized  and  plotted  on  the  theore- 
tical curve  calculated  using  the  selected  values  for  the  parameters. 
Empirically,  it  was  found  that  the  theoretical  curve  was  in  excellent 
agreement  with  the  experimental  data  points  if  the  following  values 
for  the  parameters  were  selected: 

T = 85  psec 
L/R  = 80  psec 
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Figure  HIO.  Test  configuration  for  short  circuit  current  measurement. 


and 


vU 


[ ^ /5]  i _L  = . 1 653  amperes 

I-  /r 

As  can  be  seen  from  the  plot  in  Figure  Hll,  the  time  response  of 
Eq  H20  is  in  excellent  agreement  with  the  experimental  data. 

To  evaluate  the  agreement  of  the  selected  parameters  with  the 
measured  parameters,  the  value  of  each  parameter  was  calculated  from 
the  measured  values: 

i = 2.99  m 
Q = 4 X 10"^  C 
a = .0133  m 
b = .0167  m 

p = 4tt  X 10”^  henrys/m 
0 = 2.53  X 10^  mhos/m 
L = 1 .5  yhenrys 
R = .0182  n 

The  parameters  were  calculated  to  be: 


T = = 93  ^33, 


L/R  = 82.5  psec 


r .^Q /M]  i _L 

2TT^b  ^ ^ ^ 


.1973  amperes 


These  calculated  values  are  in  agreement  with  the  theoretical 
values,  considering  the  accuracy  of  the  measurement  of  the  values 
used  in  the  calculation. 

Experimental  data  were  taken  with  additional  resistance  added  to 
the  system,  thus  decreasing  the  value  of  the  L/R  time  constant.  The 
resistance  was  increased  to  .029  ohms  and  the  theoretical  expression 
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ONE  INCH  ALUMINUM  CONDUIT  CURRENT  RESPONSE 


TIME  (MILLISECONDS) 


(Eq  H7)  was  plotted  to  fit  the  data  points;  the  value  of  L/R  was 
selected  to  be  65  psec  and  T was  selected  to  be  85  usee  to  achieve  the 
best  fit  of  the  theoretical  curve  to  the  data.  Figure  HI 2 shows  the 
plot.  The  value  of  the  inductance  was  measured  to  be  1.8  pH,  so  the 
experimental  value  of  L/R  was  62  psec.  The  theoretical  expression  is 
thus  in  excellent  agreement  with  the  experimental  data. 

The  resistance  was  then  increased  to  .11  ohms,  and  the  theoretical 
expression  was  plotted  to  fit  the  data.  The  parameters  L/R  and  T were 
selected  to  be  17  psec  and  85  psec,  respectively.  The  plot  of  the  curve 
and  the  data  is  shown  in  Figure  HI  3.  The  inductance  was  measured  to 
be  1.7  pH,  so  the  experimental  value  of  L/R  is  15.5  psec.  Again  the 
theoretical  expression  is  in  excellent  agreement  with  the  experimental 
data. 


To  demonstrate  the  effect  of  changing  the  value  of  L/R,  the  theore- 
tical curve  for  various  values  of  L/R  was  plotted  (Figure  H14). 

Summary 

The  calculated  induced  current  due  to  the  impulse  current  has  a 
relatively  simple  series  solution  which  involves  parameters  dependent 
on  conduit  dimensions  and  the  resistance  and  inductance  of  the  conduit 
sense  wire  system.  It  was  demonstrated  that  experimental  data  could 
be  put  in  the  form  of  the  theoretical  calculations  with  excellent 
agreement  for  nominal  1-in.  (25-mm)  aluminum  conduit. 
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ONE  INCH  ALUMINUM  CONDUIT  CURRENT  RESPONSE 


Figure  H12.  Theoretical  vs.  experimental  data  for  resistance  of  .029  ohms 


ONE  INCH  ALUMINUM  CONDUIT  CURRENT  RESPONSE 
T=:  85  USECS,  L/R=  17  USECS 


Figure  H13.  Theoretical  curve  vs.  experimental  data  for  resistance  of  0.11  ohms 


ONE  INCH  ALUMINUM  CONDUIT  CURRENT  RESPONSE 


TIME  (MILLISECONDS) 

Figure  H14.  Effects  of  varying  the  value  of  L/R 
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Wilson,  Monti  R.,  Transient  Current  Estimates  for  Finite  Length 
Surface  Cables,  PEM-34  (Lawrence  Livermore  Laboratory,  February  1975). 

This  report  analyzes  the  problem  of  a cable  lying  at  the  earth-air 
interface.  The  transient  response  to  a double-exponential  incident 
pulse  with  10-nsec  rise  time  is  obtained  for  various  incident  polari- 
zation, cable  lengths,  soil  conductivities,  and  terminal  impedances. 


Signal  Coupling  to  Internal  Wire 

Madle,  P.  J.,  Cable  and  Connector  Shielding  Attenuation  and  Transfer 
Irrqiedance  Measurements  Using  Quadtraxial  and  Quintaxial  Test  Methods, 

PEM-45  (Lawrence  Livermore  Laboratory,  November  1975).  This  report  des- 
cribes a frequency-dependent  transfer  impedance  and  admittance  between 
cable  shields  and  internal  conductors  which  could  also  be  applied  to 
conduit-shielded  cables.  The  following  test  fixtures  and  related  test 


procedures  are  listed:  open  wire,  coaxial  and  inverted  triaxial,  tri- 
axial,  quadraxial,  and  quintaxial. 


Vance,  E.  F.,  DNA  EMP  (Eteatvomagnetia  Pulse)  Handbook  (U),  DNA- 
2114H-2,  (Defense  Nuclear  Agency,  December  1974),  Chapter  11,  "EMP 
Coupling  to  Cables,"  ADB001204.  This  chapter  presents  a detailed  dis- 
cussion of  EMP  interaction  with  soil  and  propagation  through  the  soil. 

An  analysis  to  determine  the  induced  current  on  a base  cable  trans- 
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developed  to  guide  the  communication  system  designer  in  providing  EMP- 
resistant  facilities.  The  use  of  current-diversion  and  voltage- 
limiting  techniques  and  the  role  of  building  shields  and  facility 
ground  systems  in  EMP-resi stant  design  are  described. 

Wells,  Wm.  C.,  Nuclear  Electromagnetic  Pulse  (HEMP)  Hardened  Cables, 

Report  No.  ECOM  0341-1/ADA006642  (U.  S.  Army  Electronics  Command, 

January  1975).  Surface  transfer  impedance  Z(t)  was  determined  for  12 
prototype  cables  constructed  to  measure  the  extent  of  reduction  of  Z(t) 
possible.  Mechanisms  for  inducing  current,  antenna  effects,  and  trans- 
mission line  effects  are  discussed. 


Center,  J.  M.,  et  al..  Single-line  Modeling  of  Internal  Coupling 
for  Complex  Cable  Systems,  AFWL-TR-75-164/ADB012188  (Air  Force  Weapons 
Laboratory,  May  1976).  This  report  presents  an  analysis  of  complex 
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cables.  The  EOPRM  (Geometric  Cable  Parameter  Determination)  computer 
program  developed  for  this  application  is  described.  The  analysis  was 
conducted  for  aeronautical  systems,  but  is  also  applicable  to  fixed 
facilities. 

Cross-Coupling  in  the  SAFEGUARD  System,  Technical  Memorandum  TM- 
81  (Boeing  Aerospace  Company  for  Huntsville  Division,  U.  S.  Army  Corps 
of  Engineers,  December  1974).  Cross-coupling  models  for  small  enclo- 
sures and  conduit  runs  based  on  comparisons  with  available  test  data 
were  constructed.  These  models  were  used  to  calculate  a variety  of 
cross-coupling  situations  for  the  SAFEGUARD  system. 

Vincent,  M.  L. , C.  C.  Sutter,  and  G.  L.  Maxam,  Guideline  for  Ap- 
plication of  Single-Line  Modeling  to  B-l  Internal  Coupling,  AFWL-TR- 
75-1 35/ADBOl  1208  (A.ir  Force  Weapons  Laboratory,  April  1976).  This 
report  describes  the  development  of  the  SNGLIN  program,  in  which  coup- 
ling between  conductors  is  determined  from  the  basic  transmission  line 
equations  and  transfer  functions  between  the  conductors. 


Saturation  Effects 

Dahlen,  G. , K.  Daxberg,  L.  Hoglund,  B.  Sjoholm,  and  M.  Wik,  "A 
Survey  of  Swedish  Nuclear  Electromagnetics  (EMP).  Research,"  Intema- 
^ tional  EMC  Symposium  Record  1973,  IEEE  73  CHO  751-8  EMC  (1973),  pp 

12-18.  This  report  surveys  some  of  the  completed  and  current  research 
at  the  Research  Institute  of  Swedish  National  Defense  (FOA)  on  the 
EMP,  its  effects,  and  protection  problems.  Studies  include  tests  of 
conduits  and  EMP-shielded  cables.  A curve  showing  pulse  transmission 
through  a saturated  Fe  cable  shield  is  given. 

Kozakoff,  P.  J.,  Diffusion  of  Transient  Electromagnetic  Fields 
Thru  Saturated  Ferromagnetic  Media,  AD721906  (U.  S.  Army  Corps  of 
Engineers,  Huntsville  Division,  June  1970).  This  report  describes  a 
computer  solution  of  the  ferromagnetic  saturation  effects.  The 
solution  incorporates  material  characteristics  directly  into  the  field 
equations. 


Transmission  Line  Effects 

Brown,  Glen  L.,  General  Solutions  and  Simplified  Famulus  for 
Calculating  EMP  Induced  Currents  in  Lossless  Transmission  Lines,  PEM- 
27  (Lawrence  Livermore  Laboratory,  undated).  This  note  presents  the 
transmission  line  equation  for  EMP  coupling,  the  general  solution  for 
currents  in  the  frequency  domain,  useful  formulas  for  calculation  of 
currents  in  the  lossless  line,  and  some  experimental  comparisons.  Ex- 
amples shown  are  a transmission  line  long  in  relationship  to  pulse 
length  and  a transmission  line  short  in  relationship  to  pulse  length. 
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